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1  3,  ABSTRACT  (Maamm  200  TOCS) 

This  I'inal  report  describes  the  design  and  hardware  implementation  of  an  acousto-optic 
(AO)  multichannel  adaptive  optical  processor  (MADOP)  for  application  to  the 
c.incel lation  of  multipath  jamming  interference  in  advanced  surveillance  radars.  The 
primarv  objective  of  this  effort  was  the  design  and  preliminary  fabrication  of  a 
mu  1 1  iclianne 1  adaptive  system  that  can  perform  cancellation  of  multiple  wideband  (10  MHz] 
jamming  interference  in  the  presence  of  multipath.  The  MADOP  system  is  composed  ot 
three  primary  subsystems.  A  multichannel  t ime- integral ing  AO  correlator  performs  a 
correlation  between  the  residual  cancellation  error  and  each  of  the  auxiliary 
omnidirectional  antennas  to  arrive  at  updates  for  each  of  the  adaptive  weight  functions 
A  digital  interface  accepts  this  update  information  and  generates  the  appropriate 
.idaptive  weight  functions  for  performing  auxiliary  channel  filtering.  Finally,  a 
multichannel  AO  tapped  delay  line  (AOTDL)  filter  system  accepts  these  weight  functions  : 
through  an  AO  spatial  light  modulator  (AOSLM)  and  taps  the  auxiliary  channel  inputs  to 
form  the  estimate  of  the  noise  signal  in  the  main  receiver  channel.  This  noise  sign.il 
estimate  is  then  subtracted  from  the  main  receiver  channel  at  the  system  intermediate 
frequency  to  form  the  residual  cancellation  error  for  input  to  the  multichannel  time- 
integrating  AO  correlator,  thereby  closing  the  adaptive  loop. 
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ABSTRACT 


This  final  report  describes  the  design  and  hardware  implementation  of  an  acousto-optic  (AO) 
multichannel  adaptive  optical  processor  (MADOP)  for  application  to  the  cancellation  of  multipath  jamming 
interference  in  advanced  surveillance  radars.  The  primary  objective  of  this  effort  was  the  design  and 
preliminary  fabrication  of  a  multichannel  adaptive  system  that  can  perform  cancellation  of  multiple  wide¬ 
band  (to  MHz)  jamming  interference  in  the  presence  of  multipath.  The  MADOP  system  is  corrposed  of 
three  primary  subsystems.  A  multichannel  time-integrating  AO  correlator  performs  a  correlation  between 
the  residual  cancellation  error  and  each  of  the  auxiliary  omnidirectional  antennas  to  arrive  at  updates  for 
each  of  the  adaptive  weight  functions  A  digital  interface  accepts  this  update  information  and  generates  the 
appropriate  adaptive  weight  functions  for  performing  auxiliary  channel  filtering.  Finally,  a  multichannel  AO 
tapped  delay  line  (AOTDL)  filter  system  accepts  these  weight  functions  through  an  AO  spatial  light 
modulator  ( AOSLM)  and  taps  the  auxiliary  channel  inputs  to  form  the  estimate  of  the  noise  signal  in  the  main 
receiver  channel.  This  noise  signal  estimate  is  then  subtracted  from  the  main  receiver  channel  at  the 
system  intermediate  frequency  to  form  the  residual  cancellation  error  for  input  to  the  multichannel  time- 
integrating  AO  correlator,  thereby  closing  the  adaptive  loop. 
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1.  INTRODUCTION 


This  final  report  describes  the  design  and  hardware  implementation  of  an  acousto-optic  (AO) 
multichannel  adaptive  optical  processor  (MADOP),  for  application  to  the  cancellation  of  multipath  jamming 
interference  in  advanced  surveillance  radars.  This  Expert  in  Science  and  Engineering  (ESE)  effort  is  a 
continuation  of  an  ongoing  program  within  the  Rome  Laboratory  Photonics  Center  (References  1 
through  4).  Greater  than  half  of  the  Dynetics,  Inc.,  effort  was  performed  on-site  at  Rome  Laboratory,  Griffiss 
Air  Force  Base,  N.Y.  This  work  was  accomplished  in  conjunction  with  Photonics  Center  personnel:  CPT 
Michael  Ward,  CPT  Christopher  Keefer,  and  ILT  Harold  Andrews  under  the  In-House  Project  4600P1063; 
and  Dr  Stephen  Welstead  of  COLSA,  Inc.,  who  was  under  a  separate  ESE  contract.  All  work  reported 
herein  was  performed  under  ESE  Contract  F30602-91-D-0001 .  Subcontract  C-140545,  during  the  period  1 
February  1991  to  31  January  1992 

The  primary  objective  of  this  effort  was  the  design  and  preliminary  fabrication  of  a  multichannel 
adaptive  system  that  can  perform  cancellation  of  multiple  wideband  (10  MHz)  jamming  signals  in  the 
presence  of  multipath.  In  developing  this  system,  it  was  recognized  that  emphasis  must  be  placed  on  the 
effective  interface  of  the  MADOP  to  experimental  radar  test-beds.  For  the  radar  system  applications 
pursued  in  this  effort,  cancellation  is  best  performed  at  the  radar  intermediate  frequency  (IF),  which,  for  our 
design,  is  80  MHz.  Therefore,  the  MADOP,  when  properly  interfaced,  acts  as  a  "black  box"  component 
within  the  IF  chain  of  the  radar.  This  approach  drove  the  design  away  from  earlier  architectures  that 
performed  cancellation  at  baseband 

The  MADOP  system  is  composed  of  three  primary  subsystems:  a  multichannel  time-integrating 
AO  correlator  that  performs  a  correlation  between  the  residual  cancellation  error  and  each  of  the  auxiliary 
omnidirectional  antennas  to  arrive  at  updates  for  each  of  the  adaptive  weight  functions;  a  digital  interface 
that  accepts  the  update  information  and  generates  the  appropriate  adaptive  weight  functions  for  performing 
auxiliary  channel  filtering;  and  a  multichannel  AO  tapped  delay  line  (AOTDL)  system  that  accepts  the  weight 
functions  through  an  AO  spatial  light  modulator  (AOSLM)  and  taps  the  auxiliary  channel  inputs  to  form  the 
estimate  of  the  rvDise  signal  in  the  main  receiver  channel.  This  noise  signal  estimate  is  then  subtracted  from 
the  main  receiver  channel  at  the  system  IF  to  form  the  residual  cancellation  error  for  input  to  the 
multichannel  time-integrating  AO  correlator,  thereby  closing  the  adaptive  loop.  The  enphasis  of  the  system 
design  during  this  effort  has  been  on  the  two  AO  subsystems.  Some  assistance  has  been  provided  by 
Dynetics,  Inc.,  toward  the  design  of  the  digital  interface,  and  it  is  anticipated  that  this  digital  subsystem  will 
be  increasingly  emphasized  as  integration  into  radar  test  beds  draws  near. 

This  technical  report  is  organized  as  follows:  After  this  introduction,  the  radar  system  application 
and  design-to  requirements  are  provided  in  Section  2,  with  emphasis  on  the  integration  of  the  MADOP  into 
a  radar  testbed  at  an  IF  of  30  MHz.  This  is  followed  in  Section  3  by  a  description  of  the  adaptive  signal 
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processing  algorithms  being  Implemented,  including  initial  simulation  results.  The  analysis,  design,  and 
hardware  implementation  of  the  multichannel  AO  subsystems  constitutes  the  majority  of  the  report  and  is 
presented  in  detail  in  Section  4  Conclusions  and  rf^commendations  for  continued  development  of  this 
promising  technology  are  presented  in  Section  5 


2.  SYSTEM  CONCEPT  AND  REQUIREMENTS 


The  MADOP  has  been  designed  to  achieve  time-domain  cancellation  of  interference  sources  in 
the  presence  of  a  multipath  environment.  As  designed,  the  system  will  be  able  to  cancel  four  independent 
interference  sources,  each  having  up  to  four  multipaths,  with  each  source  having  up  to  a  1 0  MHz  bandwidth. 
The  initial  implementation  emphasizes  two  channels.  The  key  advantage  of  our  optical  approach  over  the 
current  state-of-the-art  in  digital  and  analog  electronic  implementations  is  the  ability  to  cancel  wideband 
signals  (greater  than  1  MHz)  by  using  multiple  tap  weights  per  channel.  Subsection  2.1  states  the  specific 
design-to  requirements  and  summarizes  the  radar  problem  being  addressed.  Tnis  discussion  emphasizes 
the  ultimate  integration  of  the  MADOP  black  box  into  radar  systems,  and  the  interfaces  requited. 
Subsection  2.2  reviews  the  analysis  performed  under  the  current  effort  to  assess  the  suitability  of  available 
test-beds  at  Rome  Laboratory  for  testing  the  MADOP  system. 

2.1  RADAR  SYSTEM  OVERVIEW  AND  REQUIREMENTS 

Surveillance  in  the  presence  of  multipath  jamming  interference  is  an  extremely  demanding 
problem  due  to  the  multidimensional  nature  of  the  cancellation  required.  This  is  especially  true  for  wideband 
interference  sources  that  have  several  multipath  delays.  This  subsection  will  begin  with  a  general  overview 
of  the  radar  system  block  diagram  and  how  the  MADOP  black  box  would  fit  into  the  system. 

Figure  2-1  presents  an  illustrative  summary  of  the  scenario  being  addressed.  In  this  figure,  two 
independent  interference  sources,  ni(t)  and  n2(t),  having  spectral  content  in  the  radar  system  passband, 
contaminate  the  directional  main  antenna  target  return.  s(t),  through  the  sidelobes.  The  multipath  versions 
of  the  two  interference  sources  also  contribute  a  significant  noise  term  to  the  main  channel  signal  resulting 
in  the  reception  of: 

N  M 

d(l)  =S(t)+  X  X  3201^^2 (2-1) 

n  =  1  m  =  1 

where  a^n  and  a2m  represent  relative  attenuations  due  to  multipath  losses  together  with  antenna  sidelobe 
gain  relative  to  mainlobe  gain,  and  tm  and  X2m  represent  the  signal  delays.  The  direct-path  interference  is 
given  for  n  =  1  and  m  =  1 .  It  is  assumed  that  the  interference  noise  is  on  the  order  of.  or  much  larger  than, 
the  target  return,  s(t),  resulting  in  negative  signal-to-noise  ratio  (SNR)  in  the  main  antenna.  The  two 
omnidirectional  auxiliary  antennas  receive  the  interference  noise  but  the  target  return  entering  these 
antennas  is  negligible  because  of  the  low  main-channel  SNR  assumption.  Thus,  the  two  auxiliary  antennas 
receive  the  signals: 
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2N  2M 


n  =  N  +  1 


m  =  M  +  1 


3N  3M 

N2(<)  =  Z  3inn^(t-T^^)+  ^  a2m"2('-''2m) 


n  =  2N+  1 


m  =  2M+ 1 


(2-2) 


In  addition  to  these  correlated  noises,  uncorrelated  receiver  noise  is  present  at  each  antenna. 


Figure  2-2  shows  a  simplified  block  diagram  of  a  typical  surveillance  radar  that  may  benefit  from 
the  MADOP.  This  diagram  is  not  intended  to  represent  a  specific  radar  system.  The  MADOP  is  shown  in 
this  figure  at  the  appropriate  point  in  the  signal  processing  chain.  First  note  that  the  processing  is  done  at 
the  radar  IF  after  downconversion  from  the  radio  frequency  (RF)  carrier  (3  or  10  GHz  for  example).  The 
waveform  generator  provides  the  signal  to  be  transmitted  and  also  provides  the  matched  filter  to  the  signal 
processor  to  allow  for  pulse  compression  and  pulse  integration/Doppler  filtering.  The  antenna  system, 
whether  phased-array  or  dish,  forms  a  beam  (or  multiple  independent  beams,  as  possible  with  a  phased- 
array  antenna)  that  illuminates  the  target  in  the  presence  of  interference.  The  target  reflection, 
contaminated  by  interference  and  noise,  is  received  in  the  main  channel,  while  the  interference  is  also 
received  in  auxiliary  antenna  channels  (which  can  be  portions  of  a  single  phased  array).  All  signals  are 
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downconverted  to  the  IF  and  input  into  the  MADOP.  The  output  of  the  MADOP  is  the  desired  signal  with 
the  interference  suppressed.  This  signal  then  enters  the  coherent  pulse  compression  and  Doppler  filtering 
system,  which  is  often  implemented  digitally  following  analog  to  digital  (A/D)  conversion  of  in-phase  and 
quadrature  (l&Q)  channels.  After  this  signal  processing,  target  detection  takes  place  using  some  form  of 
constant  false  alarm  rate  (CFAR)  processing.  The  data  processor  performs  such  operations  as  target 
tracking,  beam  scheduling,  system  control,  and  interface  to  the  display. 


NOTES: 

1 .  DOWN-CONVERSION  AND  A/D  CAN  BE  DONE  IN  ONE  REAL  CHANNEL 
AT  A  LOW  IF 

2  SIGNAL  PROCESSING  CAN  BE  DONE  IN  ANALOG  DOMAIN  WITH  A/D 
CONVERSION  AFTER  RANGE  AND/OR  DOPPLER  PROCESSING 


BASEBAND 


BASEBAND- IF 
MIXER 


WAVEFORM 

GENERATOR 


DOWN-CONVERSION 
IN  l&O  CHANNELS 


SIGNAL  PROCESSOR 
RANGE  COMPRESSION 
DOPPLER  PROCESSING 
DETECTION 


I 


BASEBAND 


A/D  CONVERSION 
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Figure  2-2.  Surveillance  Radar  Simplified  Schematic  with  Adaptive  IF  Processing 


As  stated  above,  the  goal  of  the  adaptive  processor  is  to  react  to  the  multipath  interference 
environnfient  in  such  a  way  that  the  main  antenna  noise  is  canceled,  resulting  in  a  satisfactory  SNR  for 
further  processing.  The  algorithm  for  achieving  this  objective,  together  with  simulated  results,  will  be 
described  in  Section  3.  The  goals  for  this  feasibility  demonstration  program,  as  coordinated  with  personnel 
in  the  Surveillance  Technology  Division  of  the  Surveillance  and  Photonics  Directorate  (OCTS),  are  shown 
in  Table  2-1.  It  is  recognized  that  although  the  loop  lock  time  (time  to  reach  steady  state)  is  critical  in  system 
applications  (500  ps  is  desirable),  the  personal  computer  interface  in  the  digital  subsystem  makes  practical 
loop  lock  times  unachievable.  The  development  of  real-time,  special-purpose  digital  interfaces  will  greatly 
increase  system  speed. 
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Table  2-1.  Feasibility  Demonstration  Goals 


System  Bandwidth  (MHz) 

10 

Number  of  Interference  Sources 

4  (initial  demo  will  accommodate  2) 

Maximum  Multipath  Delay  (ps) 

5  (corresponds  to  1 .5-km  differential  path) 

Interference  Cancellation  Ratio  (dB) 

30 

Number  of  Multipath  Delays 

>4 

Loop  Lock  Time  (ms) 

5 

Processor  IF  (MHz) 

80 

2.1.1  Signal  Representations 

The  carrier  frequency  of  the  main  and  auxiliary  antennas  is  generally  much  higher  than  the 
desired  signal  processing  IF.  Therefore,  the  radar  wilt  mix  the  signal  down  to  the  signal  processor's  IF,  f|p 
(80  MHz).  The  radar  will  not  demodulate  the  signal  down  to  baseband  and  then  modulate  the  signal  back 
to  the  IF  for  processing.  The  signal  processor  must  therefore  not  be  restricted  to  the  processing  of  double 
sideband  suppressed  carrier  (DSB-SC)  modulated  signals  (real  modulations  of  a  carrier).  An  example  of  a 
signal  that  must  be  accommodated  is  a  single-tone  interference  source.  Although  this  is  the  case,  it  is 
appropriate  to  use  OSB-SC  signals  as  test  waveforms  in  characterizing  the  processor.  In  addition,  it  is 
desirable  to  use  DSB-SC  waveforms  as  opposed  to  amplitude  modulation  (AM),  which  has  a  dc  component 
at  baseband.  These  DSB-SC  waveforms  are  more  representative  of  realistic  signals. 

In  general,  the  signal  return  and  the  interference  noise  terms  are  complex  modulations  of  a  carrier 
and  must  be  expressed  in  l&Q  channels  at  baseband.  Let  the  baseband  representation,  gn^(t),  of  a  signal, 
g(t),  be  given  as: 

g^(t)  =  g,(t)+jgQ(t)  =  A{t)exp[j4>(t)]  (2-3) 

where  gi(t)  and  gQ(t)  are  the  l&Q  channel  signals,  respectively,  and  A(t)  and  (l>(t)  are  the  signal  magnitude 
and  phase,  respectively.  On  a  carrier  frequency,  f^,  the  resultant  signal,  g(t),  is  given  as; 

g(t)  =  Re  {g^(t)exp(j27tfj.t] }  =  g,(t)cos2nfj.t-gQ(t)  sin2nfpt 

=  A(t)cos  [2ny+<|>(t)] .  (2-4) 

A  property  of  this  modulation  is  that  to  presen/e  the  information  at  baseband,  l&Q  channels  must  be 
employed.  Also,  the  bandwidth  of  a  signal  on  a  carrier  is  the  same  as  at  baseband,  where  the  I  and  Q 
channels  each  contain  half  the  bandwidth  of  the  total  signal.  For  DSB-SC  modulation,  the  carrier  modulated 
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signal  covers  twice  the  bandwidth  of  the  baseband  representation.  This  is  another  way  of  visualizing  the 
loss  of  information  when  demodulating  a  complex  signal  in  one  real  channel.  Therefore,  earlier  approaches 
to  cancellation  that  relied  on  DSB-SC  modulation  (Reference  1)  are  not  general  enough  to  handle  arbitrary 
modulations. 

To  visualize  the  correlation  that  will  result  between  the  residual  error  and  the  auxiliary  channel 
inputs  (the  algorithm  is  further  discussed  in  Section  3),  it  is  helpful  to  examine  the  bandwidth  properties  of 
signals.  The  autocorrelation  of  a  signal  will  have  a  half-power  (3  dB)  width  approximately  equal  to  the 
inverse  bandwidth.  Thus,  a  10-MHz  bandwidth  signal  will  have  an  autocorrelation  width  of  100  ns,  while  a 
1  -MHz  bandwidth  signal  will  have  an  autocorrelation  width  of  1 000  ns.  It  is  evident  that  for  multipath  delays 
of  less  than,  or  on  the  order  of  the  inverse  signal  bandwidth,  there  will  be  significant  overlap  in  the  correlation 
peaks,  which  severely  complicates  any  peak-finding  algorithm.  The  phase  of  each  correlation  may  provide 
information  useful  to  resolving  the  presence  of  two  or  more  peaks,  but  the  problem  is  still  severe.  In 
addition,  the  interaction  of  the  phases  of  the  two  overlapping  peaks  changes  the  general  appearance  of  the 
envelope  significantly.  For  example,  two  closely  spaced  peaks  can  be  made  to  look  like  one  peak,  or  a  null 
of  tens  of  dB  can  be  placed  between  them,  depending  on  the  relative  phases  of  the  peaks. 

Another  consideration  is  the  positioning  of  the  correlation  peaks  and  the  total  number  of  peaks 
for  a  given  scenario.  First,  note  that  for  a  source  and  its  multipath  replica,  the  cross-correlation  between 
two  antennas  yields  four  peaks.  The  multipath  delay  is,  in  general,  different  for  each  antenna:  if  not,  there 
appears  to  be  just  one  signal  to  cancel  without  multipath.  The  number  of  peaks  will  be  the  square  of  the 
number  of  multipaths.  The  ultimate  goal  of  the  processor  is  to  cancel  four  interference  sources,  each  having 
four  multipaths.  When  four  independent  interference  sources  are  considered,  each  having  four  multipaths, 
and  four  auxiliary  antennas,  the  total  number  of  correlation  peaks  will  be  64  per  channel  ( 1 6  per  interference 
source),  for  a  total  of  256  correlation  peaks.  These  issues  must  be  considered  in  designing  a  robust 
algorithm  for  noise  cancellation. 

2.1.2  Delay  Resolution  Requirements 

Current  program  goals  include  a  total  multipath  delay  capability  of  5  ps,  but  consideration  of  direct 
path  delay  is  also  an  important  issue.  The  issues  related  to  direct  path  cancellation  will  be  described  in  this 
subsection,  and  their  impact  on  system  development  will  be  assessed.  In  addition  to  the  impact  on  system 
design  and  performance,  the  spatial  nature  of  the  cancellation,  i.e.,  the  sidelobe  cancellation  viewpoint,  will 
become  evident. 

Figure  2-3  shows  a  scenario  containing  only  direct  path  jamming  interference  from  a  single 
source,  and  a  single  auxiliary  antenna  to  cancel  interference  in  the  main  channel.  Consider  an  antenna 
separation,  d,  and  assume  that  the  jamming  source  is  in  the  far-field.  For  this  case,  the  time-difference-of- 
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arrival  (TDOA)  of  an  interference  signal  at  the  two  antennas  is  determined  only  by  the  angle.  0,  of  the 
source  in  the  far-field.  relative  to  the  antenna  boresight.  This  TDOA  is  given  geometrically  as: 


TDOA(0)  = 


dsin0 


(2-5) 


where  c  is  the  speed  of  light.  This  forms  the  basis  for  true  time-delay  beamforming.  Figure  2-4  graphs  this 
relationship  as  the  ratio  of  the  TDOA  to  d,  the  antenna  separation.  For  example,  for  an  antenna  separation 
of  10  m  and  an  angle  off  boresight  of  5°,  the  TDOA  is  2.9  ns.  Thus,  if  the  antenna  system  is  to  cancel  the 
interference  in  this  direction  (sidelobe  cancellation  in  the  space  domain),  the  delay  must  be  set  to  precisely 
this  value.  For  two  such  interference  sources,  the  resolution  between  taps  must  be  this  precise  to  achieve 
cancellation  of  both  sources  while,  at  the  same  time,  preserving  the  mainbeam  integrity. 


Figure  2*5  further  clarifies  this  relationship  for  a  canceller  having  a  resolution  of  5  ns  over  a  5  ps 
delay  window  (1000  resolvable  spots  for  the  AOSLM  to  be  discussed  in  Section  4).  For  this  constraint,  the 
angular  resolution  of  the  canceller  will  be  a  function  of  the  antenna  separation.  Figure  2-S(a)  illustrates  the 
poor  angular  resolution  achieved  for  a  5-ns  delay  resolution  when  closely  spaced  antennas  are  employed. 
To  achieve  an  angular  resolution  of  3°,  an  antenna  separation  of  28.7  m  is  required,  as  shown  in 
Figure  2-5(b).  For  small  antenna  separations,  characteristic  of  auxiliary  antennas  that  are  part  of  a  larger 
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phased-array  main  antenna  or  mounted  on  the  main  antenna,  the  delay  resolution  must  be  very  precise  to 
cancel  direct-path  jamming  with  high  angular  resolution. 


TR  92-SI01-0002 


Figure  2-4.  TDOA  Versus  Angle  Off  Boresight  Referenced  to 
_ 1  m  of  Antenna  Separation 

if  the  interference  reaches  the  auxiliary  antenna  and  the  main  antenna  with  a  TDOA  of  less  than 
d/c  (0  =  90°),  the  appropriate  spatial  direction  will  be  cancelled.  For  this  case,  which  will  most  likely  occur 
for  each  multipath  reflection  reaching  both  antennas,  the  interference  source  appears  to  be  at  a  given 
spatial  direction.  If  the  multipath  delay  is  greater  than  this,  it  is  appropriate  to  consider  the  cancellation 
process  in  the  temporal  domain  without  consideration  of  spatial  cancellation.  In  the  special  case  of 
multipath  reaching  the  main  and  auxiliary  antennas  simultaneously,  the  system  will  have  the  undesirable 
effect  of  cancelling  the  mainbeam  return,  even  though  the  jammer  is  spatially  in  the  sidelobes. 

For  the  AOTDL  filter  system  to  be  described  in  detail  in  Subsection  4.2,  the  impact  of  this 
requirement  is  best  viewed  in  terms  of  the  number  of  resolvable  spots  in  the  AOSLM  system.  The  accuracy 
of  the  taps  can  be  very  precise,  but  the  resolution  of  two  closely  spaced  taps  will  most  likely  be  limited  to 
less  than  the  time  bandwidth  (BT)  product  of  the  AOSLM,  typically  on  the  order  of  1000.  A  possible  solution 
to  this  constraint  is  to  use  very  high  temporal  resolution  over  a  very  small  delay  window,  and  then  to  form 
coarser  taps  over  the  total  5  ps  window.  This  coukJ  possibly  be  achieved  in  a  single  optical  system,  or  may 
be  performed  by  a  large  delay  window  multipath  optical  canceller  in  conjunction  with  a  traditional  electronic 
skjelobe  canceller. 
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ANGULAR  RESOLUTION  (deg)  ANGULAR  RESOLUTION  (deg) 


MAIN-AUXILIARY  ANTENNA  SEPARATION  (m) 
(a)  0  to  20  m 


MAIN-AUXILIARY  ANTENNA  SEPARATION  (m) 
(b)  20  to  200  m 

TR-92-SI01  0003 


Figure  2-5.  Angular  Resolution  Versus  Antenna  Separation 
Referenced  to  5-ns  Delay  Resolution 
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SURVEILLANCE  TEST-BED  ASSESSMENT 


During  the  on-site  support,  meetings  with  personnel  from  organizations  outside  the  Photonics 
Center  were  held  to  identify  potential  test-beds  for  the  MADOP  system.  The  first  meeting  was  concerned 
with  the  L-band  and  S-band  radars  operated  by  OCXS.  The  discussions  were  specifically  concerned  with 
attempting  to  determine  if  the  MADOP  would  be  useful  in  the  aforementioned  radars  and  whether  one  of 
the  radars  could  serve  as  a  test-bed  for  evaluating  the  MADOP.  It  was  concluded  that  the  mission  of  these 
two  radars  is  such  that  the  MADOP  would  not  be  useful.  In  particular,  these  radars  appear  to  be  test-beds 
that  are  used  to  evaluate  concepts  proposed  by  other  organizations.  OCTS  would  not  have  a  need  to 
incorporate  an  adaptive  interference  canceller  since  it  would  not  enhance  the  role  of  the  radars  as  test-beds. 
The  L-band  and  S-band  radars  also  would  not  serve  as  good  test-beds  for  evaluating  the  MADOP.  This 
conclusion  stems  from  a  perception  that  it  would  be  very  difficult  to  interface  the  MADOP  to  the  radars  or  to 
configure  them  to  provide  data  that  could  be  used  to  exercise  the  MADOP  in  an  off-line  fashion.  The  main 
difficulty  lies  in  the  fact  that  there  was  no  evidence  that  the  radars  had  the  auxiliary  antennas  and  receiver 
channels  required  for  this  application. 

A  second  meeting  was  held  with  personnel  from  OCDR.  The  purpose  of  this  meeting  was  to 
discuss  the  C-band  radar  operated  by  OCDR.  This  radar  is  a  6  by  12  ft  planar  phased  array  that  was 
developed  for  electronic  counter-countermeasure  (ECCM)  testing.  As  with  the  L-  and  S-band  radars,  the 
mission  of  this  radar  is  not  such  that  the  MADOP  would  be  useful  in  it.  However,  it  could  serve  as  a  test¬ 
bed  for  evaluating  the  MADOP.  It  has  an  auxiliary  antenna  and  extra  receiver  channels  that  can  be  used 
to  interface  to  the  MADOP.  All  the  receiver  channels  are  readily  available  and  provide  outputs  at  the  proper 
IF.  They  also  appear  to  have  bandwidths  that  are  compatible  with  the  capabilities  of  the  MADOP.  The  radar 
has  a  second  auxiliary  antenna  that  could  be  mounted  near  the  existing  auxiliary  antenna  to  provide  two 
auxiliary  channels.  Furthermore,  it  might  be  possible  to  form  additional  auxiliary  channels  by  accessing 
individual  subarrays  of  the  main  array.  However,  before  such  an  approach  is  tried,  it  would  be  necessary 
to  obtain  details  on  the  scanning  characteristics  of  the  subarray.  If  the  subarrays  are  scanned,  this  could 
introduce  complications  into  the  operation  of  the  MADOP.  Other  attractive  features  of  the  C-band  array 
radar  are  the  two  jammers  and  a  boresight  tower  associated  with  it.  If  a  corner  reflector  or  some  type  of 
signal  repeater  could  be  attached  to  the  boresight  tower,  it  could  provide  a  target  signal  source  for  some  of 
the  testing. 

A  third  meeting  was  held  with  personnel  from  C3BA.  This  group  has  an  X-band  communications- 
type  system  that  they  use  to  study  sidelobe  cancellation  and  other  functions.  The  system  consists  of  a 
12-horn  antenna  with  one  receiver  channel  per  horn.  All  of  the  beamforming,  sidelobe  cancellation,  and 
other  processing  are  done  at  baseband  by  a  digital  computer.  They  have  not  pursued  the  area  of  noise 
cancellation,  nor  are  they  conveniently  set  up  to  do  so.  Because  of  this,  it  is  questionable  as  to  whether  this 
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system  would  provide  a  good  test-bed  for  the  MADOP.  If  the  MADOP  were  to  be  modified  to  serve  as  a 
sidelobe  canceller,  this  facility  might  be  useful  for  testing. 

A  fourth  meeting  was  held  with  personnel  from  OCDE  on  the  potential  application  of  photonics  to 
radars.  A  system  that  was  of  some  interest  in  terms  of  the  MADOP  was  their  256-element  patch  array.  They 
apparently  plan  to  use  this  array  to  study  sidelobe  cancellation  and  adaptive  beamforming  techniques.  At 
a  later  time,  this  radar  could  also  serve  as  a  test-bed  for  the  MADOP.  Another  thought  is  that  if  the  MADOP 
was  converted  to  a  sidelobe  canceller  or  a  piece  of  hardware  to  support  adaptive  beamforming,  OCDE 
might  be  interested  in  using  it  in  their  patch  array. 

It  appears  that  the  MADOP  is  maturing  to  a  point  where  it  can  be  tested  in  a  radar.  To  this  end, 
there  are  several  issues  related  to  potential  performance  and  system  implementation  that  must  be 
addressed.  One  of  these  is  a  concern  for  how  the  MADOP  performs  relative  to  a  multichannel  jamming 
suppression  system  that  does  not  contain  the  delay-line  feature.  A  further  concern  is  whether  or  not  the 
existing  version  of  the  MADOP  offers  better  performance  than  an  interference  canceller  that  does  not 
include  the  delay-line  feature.  There  is  a  related  issue  of  whether  the  nondelay-line  version  could  be  easily 
implemented  with  optical  processing  techniques.  It  is  recommended  that  this  area  be  pursued  with  analysis 
and  simulation. 

Another  area  of  performance  that  should  be  considered  is  the  convergence  time  of  the  MADOP. 
This  impacts  the  implementation,  use  and  performance  of  the  radar  when  it  incorporates  the  MADOP.  It 
may  be  that  all  pulses  transmitted  during  the  convergence  period  will  not  be  available  for  other  signal 
processing  functions.  The  net  effect  of  this  is  a  loss  in  radar  resources.  If  the  convergence  period  is  a 
significant  portion  of  the  processing  interval,  this  could  seriously  impact  system  performance.  Another  area 
related  to  convergence  time  is  the  use  of  the  MADOP  in  a  continuously  scanning  search  radar.  If  the 
convergence  time  is  too  large,  the  scan  rate  of  the  radar  would  have  to  be  slowed,  which  could  affect  overall 
performance  of  the  radar.  Of  course,  if  jamming  is  severe  enough,  the  resultant  performance  degradation 
could  be  much  worse  due  to  convergence.  Also  related  are  the  effects  of  MADOP  operation  on  clutter 
cancellation,  pulse  distortion,  target  parameter  measurements,  and  probably  many  others. 

In  the  area  of  system  implementation,  the  problem  exists  as  to  just  how  the  jammer  suppression 
system  should  be  irtipiemciited  in  the  sum  and  difference  channels  of  the  radar  and  what  effect  the  jammer 
suppression  system  would  have  on  monopulse  performance.  It  should  also  be  noted  that  for  search  radars, 
which  typically  do  not  employ  monopulse  techniques,  this  may  not  be  a  problem;  for  track  radars  however, 
it  will  definitely  be  an  issue.  Furthermore,  the  impact  on  performance  will  depend  upon  the  type  of 
monopulse  processing;  e  g.,  two  channel,  three  channel,  or  four  channel.  There  is  also  a  need  to  consider 
the  implementation  impact  on  multiple-beam  antennas. 


2-10 


3.  ADAPTIVE  CANCELLATION  ALGORITHM  AND  SIMULATION 


The  cancellation  of  noise  from  desired  signals  using  auxiliary  channel  inputs  has  been  applied  to 
a  number  of  systems  with  success  (References  5  through  7).  In  addition,  the  application  of  optical,  and 
more  specifically  acousto-optic  (AO),  signal  processing  techniques  to  the  adaptive  filtering  problem  has 
been  researched  by  others  (References  8  through  14).  Subsection  3.1  will  present  a  summary  of  the 
specific  adaptive  algorithm  being  implemented  by  the  MAOOP  system,  and  will  point  out  differences 
between  this  algorithm  and  those  employed  by  others.  Following  this  theoretical  development,  a  number  of 
alternative  weight  update  and  application  schemes  will  be  compared  through  simulation  in  Subsection  3.2. 

3.1  ADAPTIVE  CANCELLATION  ALGORITHM 

The  standard  multidimensional  tapped  delay  line  implementation  of  the  least-mean-square  (LMS) 
algorithm  is  shown  in  Figure  3-1.  We  note  that  the  LMS  algorithm  is  based  on  steepest  descent  techniques, 
where  the  weight  function  is  updated  in  the  direction  opposite  the  gradient  at  the  current  weight  value.  This 
will  cause  the  weights  to  move  in  such  a  way  as  to  reduce  the  mean-squared-error  (MSE)  along  the  steepest 
gradient.  The  input  signals  from  the  auxiliary  antennas,  N,(t),  are  convolved  with  the  weight  functions  W|(t;T), 
which  are  functions  of  time,  t,  and  delay,  x,  to  yield  an  estimate,  dest(f).  of  the  main  channel  signal: 

P  2 

^  X  1  Wj(t;x)Nj(t-T)dT  (3-1) 

i  =  i 

2 

where  P  is  the  number  of  auxiliary  antennas  and  Tg  is  the  length  of  each  tapped  delay  line.  This  estimate 
is  subtracted  from  the  main  channel  signal,  d(t),  to  form  an  error  signal.  e(t).  The  weights  are  then  updated 
according  to  the  degree  of  correlation  between  this  error  and  the  auxiliary  channel  inputs.  The  i’^  weight 
function  is  represented  mathematically  as: 

t 

W|(t;T)  =  W|(0;x)  +  Ja(t)e(t)Nj*  {t-x)dt  (3-2) 

0 

where  a(t)  is  a  time-varying  (in  general)  acceleration  parameter.  As  this  process  continues  in  time,  the 
adaptively  changing  weights  converge  to  a  steady-state  solution  (if  the  signal  environment  is  stationary). 
This  cancels  the  noise  in  the  main  channel,  which  equivalently  yields  the  minimum  MSE  and  maximum 
output  SNR. 
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Figure  3-1.  Classical  Multichannel  LMS  Adaptive  Filter 


This  classical  LMS  algorithm  performs  a  weight  update  at  each  point  in  time.  If  a  discrete 
algorithm  is  employed,  the  weight  update  is  achieved  at  each  time  step,  which  must  be  at  a  frequency 
greater  than  twice  the  widest  bandwidth  to  be  processed.  Thus,  for  10  MHz  bandwidth  signals,  the  time 
step  and  weight  update  occur  every  50  ns.  Figure  3-1  implies  that  taps  must  be  positioned  every  50  ns  if 
discrete  taps  are  employed.  The  integration  shown  in  Equation  3-2  occurs  over  the  entire  observation 
period,  resulting  in  a  continuous  weight  update.  In  addition,  as  the  weights  are  updated,  the  nrror  continues 
to  be  reduced,  which  then  feeds  back  into  the  weight  calculation  integration. 

The  algorithm  being  implemented  by  the  MADOP  differs  from  this  classical  LMS  algorithm.  The 
fitter  function  given  by  Equation  3-1  is  implemented  in  the  classical  form,  but  the  weight  update  given  in 
Equation  3-2  is  implemented  differently.  For  our  system,  the  correlation  between  the  error  and  the  auxiliary 
antenna  signals  is  performed  over  some  window  of  time,  T.  during  which  the  weight  functions  do  not 
change.  This  window  of  time  has  typically  been  greater  than  2  ms  for  the  single  channel  system,  and  is 
ultimately  limited  to  approximately  1  ps  due  to  linear  detector  array  readout  times  and  limited  detector 
sensitivity.  After  this  correlation  has  been  obtained,  the  result  is  used  to  update  the  weight  functions,  a 
process  that  takes  an  additional  increment  of  time,  t^.  Thus,  Equation  3-2  is  modified  to  yield; 
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(3-3) 


M  i(T  +  U+T 

Wj[M(T  +  t^j)  :-cl  =  Wj(0;T) -t- ^  J  a.e(t)Nj*  (t-T)dt 

)  =  ■'  i(T  +  U 

where  M  is  the  processing  iteration  step.  If  t^.  the  weight  update  time,  is  not  less  than  or  equal  to  T,  then 
during  the  time  that  the  weights  are  being  updated  there  is  no  consideration  of  the  signal  environment  (no 
correlations  are  being  integrated  in  Equation  3-3).  The  incremental  changing  of  the  weights  and  the  loss  of 
correlation  information  need  to  be  further  studied  to  determine  the  impact  on  system  performance. 

This  windowed  algorithm,  as  described  in  Equation  3-3,  is  similar  to  the"  block"  adaptive  filtering 
approach,  which  has  been  increasingly  researched  in  the  digital  adaptive  filtering  community 
(References  15  through  17).  The  primary  motivation  for  this  interest  is  the  compatibility  of  block  routines 
with  fast  convolution  algorithms  that  employ  fast  Fourier  transforms  (FFT).  In  these  adaptive  filters,  the 
weights  are  updated  based  on  a  window  of  data,  rather  than  at  each  time  step.  This  is  characteristic  of  the 
MADOP  algorithm.  An  interesting  result  related  to  these  block  routines  is  that  the  convergence  generally 
proceeds  much  faster  than  for  the  classical  LMS  routine  due  to  the  reduction  in  the  gradient  estimate  noise. 
Also,  highly  correlated  input  data  (resulting  in  a  high  correlation  matrix  condition  number  and  generally  slow 
convergence)  are  better  processed  through  block  techniques. 

Another  algorithmic  consideration  is  the  need  to  develop  complex  weights  based  on  complex 
correlation  outputs.  The  emphasis  to  date  in  the  MADOP  program  has  been  the  generation  of  real  weights. 
If  the  multipath  summations  described  in  Equation  2-2  are.  in  fact,  strictly  time  delays  with  no  phase  shift 
imposed  on  the  RF  (and  IF)  carrier,  and  if  there  are  no  propagation  or  receiver-induced  phase  shifts,  then 
the  complex  weighting  scheme  will  not  be  required.  For  other  situations,  performance  may  be  degraded  if 
complex  weights  are  not  employed. 

An  additional  algorithm  is  being  considered  that  will  provide  a  discrete  number  of  taps  based  on 
the  location  of  the  correlation  peaks  arising  in  Equation  3-3  (Reference  18).  Several  issues  must  be 
addressed  relating  to  this  approach.  First,  the  number  of  peaks  arising  in  the  correlation  window  grows  very 
large  in  complex  signal  environments,  as  described  at  the  end  of  Subsection  2.1 .1 .  Also,  for  a  single  source 
with  two  multipaths,  the  number  of  correlation  peaks  in  the  cross-correlation  will  be  four,  and  there  remains 
the  ambiguity  as  to  which  peak,  or  peaks,  to  choose.  When  single  tones  or  combinations  of  single  tones  (a 
two-tone  OSB-SC  signal  for  example)  are  interference  sources,  there  are  a  number  of  peaks  because  of  the 
periodicity  of  the  correlation.  Also,  a  signal  composed  of  the  sum  of  sinusoids  of  the  same  frequency  is  still 
a  sinusoid  at  that  same  frequency,  with  an  amplitude  and  phase  shift,  so  it  is  inappropriate  to  consider 
multipath  with  sinusoidal  inputs.  Note  that  this  comment  also  applies  to  DSB-SC  two-tone  signals  created 
by  modulating  a  sinusoid  onto  a  carrier. 
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3.2 


SIMULATION  ACTIVITIES 


Three  software  simulations  of  the  MADOP  system,  known  as  Optical  Jammer  Cancellers  1,  2, 
and  3  (OJC1 ,  OJC2,  OJC3),  were  developed  to  assess  alternative  weighting  schemes.  These  simulations, 
developed  using  MATLAB  software  obtained  from  MathWorks,  Inc.,  are  provided  in  Appendix  A  of  this 
report.  The  three  simulations  employ  a  pulse  input  signal  immersed  in  white  noise  (band-limited  in  OJC3). 
The  three  simulations  differ  in  their  weight  schemes  as  described  below. 

3.2.1  QJSl 

OJC1  demonstrates  the  capability  to  cancel  wideband  interference  after  several  iterations  using 
any  of  three  weighting  schemes,  namely  a  thresholded  weight  vector  update,  a  peak  weight  update 
selection,  and  a  full  weight  vector  update.  The  fully-commented  code  appears  in  Appendix  A,  where  the 
"%"  symbol  denotes  a  comment  in  the  software  Parameters  used  in  the  software  are  summarized  in  the 
Appendix  in  Table  A-1,  and  a  flowchart  of  the  software  for  OJC1  is  shown  in  Figure  A-l.  Figure  3-2  was 
generated  using  the  thresholded  weight  vector  update  scheme,  which  is  not  commented  out  in  Appendix  A. 

The  primary  motivation  for  simulating  the  thresholded  weight  vector  schemes  is  the  desire  to 
minimize  the  number  of  taps  applied  in  the  AOTOL,  since  this  drives  the  complexity  of  the  electronic 
interface  for  weight  vector  application.  This  is  done  in  either  of  two  ways.  The  thresholded  weigW  vector 
update  normalizes  delwxo  the  range  [0. 1]  and  sets  to  zero  any  value  greater  than  0.1  (this  value  is  easily 
changed  to  simulate  other  thresholds)  The  peak  weight  update  scheme  selects  the  largest  unnormalized 
weight  and  sets  all  other  weights  to  zero.  In  both  cases,  the  weight  update  is  added  to  the  previous  weight 
vector  to  generate  the  new  weight  vector,  which  is  not  further  thresholded. 

The  acceleration  paranieter  a  was  chosen  on  recommendations  from  Dr.  Welstead.  If  the 
previous  weight  update  vector  delwold  is  zero,  or  if  there  is  no  change  between  delwold  and  the  current 
weight  update  vector  delw,  then  a  =  0  5  if  this  is  not  the  case,  the  previous  value  of  a  {aoldi  is  changed 
using  the  equation 

.  (delwold) '(delwold)^ 

a  =  aold— - — - - — =  (3-4) 

(delwold)  '(delwold-delw) 

where  the  superscript  T  denotes  the  transpose  operation.  Thus,  the  numerator  and  denominator  are  vector 
inner  products.  This  equation  states  that  (or  fast  changing  weight  update  vectors  (large  denominator),  a 
decreased  and  vice-versa  This  helps  maintain  a  fast  rate  of  convergence  while  minimizing  the  possibility 
of  unstable  operation. 
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Figure  3-2.  OJCI  Simulation  Performance 
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Figure  3-2.  OJCl  Simulation  Performance  (Continued) 
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The  simulated  signal  environment  has  inherent  simplifications  that  will  be  upgraded  in  future 
simulation  developments.  The  spacing  between  time  samples  is  0.5  ps  (20  MHz  sampling  rate)  and  the 
simulated  signal  environment  is  periodic  with  a  period  of  10  ps  (200  samples).  The  uncorrupted  received 
signal  s  is  a  20-ps  pulse  (500  kHz  bandwidth)  and  the  jamming  signal  is  a  Gaussian-distributed  zero-mean 
random  vector  having  a  bandwidth  of  20  MHz  (equal  to  the  sampling  rate).  The  sum  of  the  received  and 
jamming  signals,  with  a  known  signal-to-jammer  power  ratio,  represents  the  received  main  channel  signal 
d.  The  auxiliary  channel  is  equivalent  to  the  jamming  signal  in  the  main  channel  but  delayed  by  1  ps  and 
weighted  by  the  auxiliary-to-main  channel  voltage  ratio  AM{AM=  1  in  OJC1  but  can  easily  be  changed).  In 
the  current  simulation  the  signal  and  jammer  representation  are  constant  during  the  adaptive  cancellation 
algorithm  and  thus  represent  a  very  idealized  situation.  This  will  be  changed  to  a  more  realistic  signal  and 
jammer  representation  in  future  studies. 

The  function  filter(b.  a,  x)  filters  a  vector  x  with  the  tapped  delay  line  filter  described  by  the 
difference  equation 

y(n]  =  b(1)x(nl  -t-b(2)x(n-1]  +  ...  +b{njj+1)x(n-njj] 

-a (2) yin- 11  -  ...-a(n3+ 1)y (n-ngj 

where  b  is  a  vector  of  length  nb  corresponding  to  the  weight  vector  in  our  analysis,  and  the  coefficients  a 
are  all  set  to  zero  except  a(1)  =  1  (the  filter  is  nonrecursive  in  our  analysis).  In  the  z-transform  frequency 
domain  this  function  is  expressed  as 

b(1)  +  b(2)z‘^  +  ...  +  b(n.  +  l)z 

Y(z)  =  - ^ - _X(z).  (3-6) 

1  +  a(2)z'^  +  ...  +a(n3+  l)z 

a 

Therefore,  in  OJCI  the  auxiliary  channel  jammer  signal  na  is  filtered  by  the  finite  impulse  response  (FIR) 
filter  having  weights  w  using  the  equation 

y  =  filter  (w,  1,  na) .  (3-7) 

The  function  xcorr(x,  y,  “biased")  forms  the  biased  estimate  of  the  cross-correlation  of  the  vectors  x  and  y 
given  by 


n  -  |m|  -  1 

C(m)  =  1  Y,  x(n)y*  (n-^m).  (3-8) 

n  =  0 

The  plots  shown  in  Figure  3-2  demonstrate  clearly  the  effectiveness  of  the  OJCl  algorithm  for 
this  signal  environment.  After  five  iterations,  the  jammer  power  in  the  main  channel  has  been  reduced  from 
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10.7  dB  to  -17.5  dB,  an  improvement  of  28.2  dB,  as  shown  in  Figure  3-2  (e).  The  pulse,  which  is  initially 
immersed  in  the  noise,  is  clearly  visible  in  the  error  output  (the  best  filtered  version  of  s).  Note  also  that  the 
acceleration  parameter  a  increases  as  the  noise  is  suppressed,  allowing  for  faster  convergence  while 
maintaining  stability. 

3.2.2  OJC2 

The  routine  OJC2  is  similar  to  OJCI  with  the  addition  of  a  multipath  replica  of  the  jammer.  The 
main  channel  receives  two  replicas  of  the  jammer  separated  by  10  samples  (0.5  ps)  in  addition  to  the  pulse 
s.  The  auxiliary  channel  contains  the  same  two  replicas  but  this  time  separated  by  25  samples  (1 .25  ps). 
Because  the  multipath  environment  is  more  severe  for  this  case,  this  single  auxiliary  channel  system  does 
not  cancel  this  multipath  jamming  signal  adequately,  as  is  clearly  evident  in  Figure  3-3.  Even  after  1 1 
iterations,  the  noise  still  dominates  the  s  signal. 

3.2.3  aica 

The  routine  OJC3  is  similar  to  OJCI  except  that  the  jamming  signal  is  now  passed  through  a 
2-MH2,  5-pole  Butterworth  filter.  Also,  for  the  results  shown  in  Figure  3-4,  the  entire  weight  update  vector 
was  employed,  as  shown  by  the  use  of  comments  in  the  OJC3  code  in  Appendix  A.  It  is  observed  once 
again  that  the  system  effectively  cancels  the  jammer  as  shown  after  five  iterations  in  Figure  3-4  e.  For  this 
bandlimited  noise  case,  the  improvement  in  jammer  power  in  the  main  channel  is  18.8  dB. 
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Figure  3-3.  OJC2  Simulation  Performance 
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Figure  3-4.  OJC3  Simulation  Performance 
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Figure  3-4.  OJC3  Simulation  Performance  (Continued) 
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4.  MULTICHANNEL  ADAPTIVE  OPTICAL  PROCESSOR  DESIGN  AND 
IMPLEMENTATION 

The  use  of  AO  signal  processing  techniques  for  adaptive  filtering  has  been  studied  by  others.  For 
example,  time-domain  implementations  employing  optical  integrators  such  as  liquid  crystal  light  values 
(LCLV)  and  photorefractive  crystals  have  been  demonstrated  (References  8,9,10,  and  11).  An  architecture 
not  requiring  an  integrating  device  was  proposed  in  Reference  12.  VanderLugt  has  developed  an 
equivalent  frequency-domain  architecture  that  also  does  not  require  an  optical  integrator  (Reference  13). 
An  excellent  review  and  performance  comparison  of  several  of  these  architectures  is  provided  in 
Reference  14. 

As  described  in  Subsection  3.1 ,  the  MADOP  architecture  implements  the  block  LMS  algorithm  as 
opposed  to  the  classical  LMS  approach.  This  is  achieved  through  the  integration  of  the  weight  update 
correlation  information  (over  a  finite  block  of  data)  on  a  charge  coupled  device  (CCD)  detector  array.  This 
result  is  then  used  to  update  the  adaptive  weight  vector  that  is  applied  to  the  AOTDL  filter  to  generate  the 
desired  main  channel  estimate.  This  architecture  is  also  being  implemented  using  multichannel  AO  cells 
and  parallel  processing  of  multiple  auxiliary  channel  inputs.  This  multichannel  processor  is  the  first  such 
architecture  reported. 

Specific  accomplishments  made  on  the  optical  implementations  of  the  time-integrating  correlators 
for  weight  update  calculation  include; 

1 .  Modification  of  an  existing  single-channel,  in-line,  time-integrating  correlator  (Reference  1 )  to 
accommodate  coherent  processing  at  the  system  IF; 

2.  Performance  characterization  of  the  modified  in-line  architecture  using  laboratory-generated 
test  waveforms; 

3.  Design  and  development  of  a  two-channel,  two-path  Mach-Zehnder  time-integrating 
correlator; 

4.  Performance  characterization  of  the  two-path  architecture  using  laboratory-generated  test 
waveforms; 

5.  Comparison  of  the  Mach-Zehnder  two-path  architecture  to  the  in-line  architecture;  and 

6.  Definition  of  post-detection  operations  for  the  interface  between  the  weight  update  calculation 
and  the  tapped  delay  line  filter,  and  verification  of  initial  post-detection  processing  routines 
developed  by  ILT  Andrews. 

Specific  accomplishments  made  on  the  optical  implementations  of  the  AO  tapped  delay  line  filter  include: 
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1 .  Design  and  mathematical  description  of  a  coherent  AOTDL  architecture  employing  a  Mach- 
Zehnder  interferometer,  a  modification  of  a  prior  baseband  AOTDL  architecture 
(Reference  1): 

2.  Fabrication  of  a  single-channel  version  of  this  system; 

3.  Test  and  analysis  of  the  coherent  architecture  that  employs  a  Mach-Zehnder  interferometer 
to  provide  an  output  signal  at  the  desired  IF  without  the  effects  of  the  frequency  shift  induced 
bythe  AOSLM; 

4.  Development  of  a  coherent  noninterferometric  architecture  with  an  electronic  reference  that 
overcomes  the  true  time  delay  shortcomings  of  the  Mach-Zehnder  architecture; 

5.  Fabrication  of  a  two-channel  version  of  this  processor  at  80  MHz  center  frequency,  and 

6.  Initial  performance  characterization  of  this  architecture  using  laboratory-generated  test 
waveforms. 

4.1  OVERALL  SYSTEM  ARCHITECTURE 

The  design  and  implementation  of  these  AO  subsystems  is  the  subject  of  this  section.  First  the 
overall  system  will  be  described,  followed  by  the  theory,  design,  and  test  results  for  the  AO  time-integrating 
correlators  and  the  AO  tapped  delay  line  filter  architectures. 

The  adaptive  IF  canceller  shown  in  Figure  2-2  accepts  the  main  channel  signal  and  auxiliary 
channel  inputs,  and  outputs  the  filtered  main  channel  signal  to  the  radar  signal  processor.  The  MADOP 
implementation  of  the  adaptive  IF  canceller  employing  two  auxiliary  channels  is  shown  in  Figure  4-1.  The 
AO  time-integrating  correlator  accepts  the  two  auxiliary-channel  inputs  and  correlates  these  against  the 
error  signal  (the  filtered  main-channel  signal).  The  two  finite-time  correlations  that  result  represent  the 
weight  update  information  and  are  input  into  a  personal  computer  (PC)  to  generate  the  adaptive  weight 
vectors.  These  two  sets  of  weights  are  then  applied  to  the  AO  tapped  delay  line  architecture  to  filter  each 
of  the  auxiliary  channel  inputs  to  arrive  at  a  multichannel  estimate  of  the  noise  in  the  main  channel 
(assuming  low  SNR  in  the  main  channel  and  therefore  no  target-return  signal  information  in  the  auxiliary 
channel  inputs).  This  noise  estimate  is  subtracted  from  the  main  channel  signal  to  form  the  error  signal  to 
be  input  into  the  time-integrating  correlator,  thereby  closing  the  loop. 

Potential  signal  and  interference  sources  for  testing  and  characterizing  the  MADOP  in  a 
laboratory  environment  are  shown  in  Figure  4-2.  The  Lecroy  9100  Arbitrary  Function  Generator  (AFG)  can 
be  used  to  generate  baseband  waveforms  such  as  tones,  linear  frequency  modulation  (LFM),  and  noise¬ 
like  modulations,  for  interference  source  simulation.  These  baseband  signals  can  then  be  mixed  up  to  the 
80-MHz  IF  using  DSB-SC  modulation.  Another  option  is  to  employ  a  noise  generator  that  outputs  a 
waveform  at  each  80  MHz  IF,  thus  requiring  no  mixing.  The  signal  source  can  be  provided  using  the  Hewlett 
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Packard  (HP)  Frequency  Agile  Signal  Simulator  (PASS)  or  the  General  Electric  (GE)  Radar  Signal  Simulator 
(RSS),  both  of  which  are  available  in  the  Photonics  Center.  Another  possibility  is  to  generate  a  third 
uncorrelated  interference  source  using  the  options  shown  in  Figure  4-2  (a)  and  call  this  third  source  the 
signal. 
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_  Figure  4-1.  MADOP  System  Layout  for  Tivo  Channels 


The  electronic  interfaces  between  the  two  AO  subsystems  are  shown  in  Figure  4-3.  Referring  to 
Figure  4-3  (a),  the  outputs  from  the  AO  time-integrating  correlator  are  centered  at  the  carrier  frequency  of 
the  correlation  information  after  detector  readout.  A  Girard  3197  interlace  board  then  digitizes  the  data  with 
8-bit  resolution  and  stores  the  digital  data  in  memory  for  further  digital  processing.  The  weights  are  updated 
in  the  PC  and  provide  the  AOTDL  filter  weights.  In  the  future,  the  PC  will  provide  automatic  gain  control 
(AGC)  to  the  error  signal  to  maintain  strong  correlation  results  as  the  error  is  adaptively  reduced. 

The  error  generation  interface  is  shown  in  Figure  4-3  (b).  The  photodiode  output  is  at  a  center 
frequency  equal  to  the  system  IF  minus  the  electronic  reference  frequency,  which  is  approximately  15  MHz 
(see  Subsection  4.3.2  for  details  on  the  electronic  reference).  After  amplification  this  signal  is  mixed  back 
up  to  the  80-MHz  IF,  bandpass  filtered,  and  subtracted  from  the  main  channel  signal.  This  error  signal  then 
drives  the  time-integrating  correlator,  which  amplifies  and,  in  the  future,  performs  AGC  for  matching  signal 
levels  in  the  correlator  system. 
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Figure  4-3.  Electronic  Interfaces 
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TIME-INTEGRATING  CORRELATION  FOR  WEIGHT  FUNCTION  CALCULATION 


As  mentioned  previously,  the  time-integrating  correlator  portion  of  the  MADOP  is  used  for  the 
calculation  of  the  weight  function  and  generation  of  the  inpufs  to  the  AOTDL  filter.  We  consider  the  error 
signal,  e(t),  and  auxiliary  antenna  input,  N(t),  given  at  baseband  as  the  corrplex  modulations 

e(t)  =  Ag(t)exp[i<D^(t)l  (4-1) 


and 


N(t)  =  A^(t)exp[i<t»^(t)l.  (4-2) 

The  coherent  (magnitude  and  phase)  correlation,  as  described  in  Equation  3-2  with  a  incorporated  in  e(t), 
is  given  by; 


T 

R  (X)  =  jAg  (t)  exp  (j<|)^  (t)  ]  A^  (t  -  X)  exp  [- j<(>^  (t  -  x)  1  dt  =  Ap  (x)  exp  [j(t)p  (t)  ]  (4-3) 

0 

where  Ar(x)  and  <t)R(x)  are  the  correlation  magnitude  and  phase,  and  T  is  the  total  integration  time.  This 
integral  can  be  realized  on  the  CCD  array  (all  values  are  positive  and  real)  by  modulating  R(x)  onto  a  spatial 
carrier,  f^,  and  adding  a  bias.  Thus,  the  desired  correlatton  output  is 


R(x)  =  bias  +  Re 


exp  (j27cf.jX]  je  (t)  N*  (t  -  X)  dt 


=  bias  Ap  (X)  cos  [2n\^i  +  <t!p  (x)  1 


(4-4) 


which  can  be  coherently  demodulated  to  obtain  the  magnitude  and  phase  of  the  correlation.  In  the  following 
subsections,  techniques  for  achieving  this  coherent  correlation  operation  in  a  multichannel  configuration  will 
be  described.  In  addition,  the  design,  testing,  and  initial  performance  evaluations  of  the  different  time- 
integi’afing  correlator  architectures  will  be  presented. 


4.2.1  Design  Alternative  I:  In-Line  Multichannel  Correlator  Architecture 


Figure  4-4  illustrates  the  conceptual  in-line  correlator  architecture  for  a  two-channel 
configuration.  The  actual  optical  layout  is  depicted  later  in  Figure  4~7(a).  A  modulated  laser  diode  is  used 
as  the  optical  source  for  this  architecture.  Also  shown  in  the  figure  is  a  diagram  of  the  spectral  inputs  to  the 
laser  diode  and  multichannel  AO  cell  for  coherent  operation.  As  shown  in  the  figure,  beam  conditioning  and 
replicating  optics  are  used  in  conjunction  with  the  laser  diode  to  illuminate  separate  channels  of  the 
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Figure  4-4.  Conceptual  Layout  for  In-Line  Correlator 


multichannel  AO  cell.  The  diffracted  signals  from  the  AO  channels  are  schlieren  imaged  onto  separate  CCD 
arrays  where  they  are  time  integrated  and  read-out  to  calculate  the  weight  updates.  Schlieren  imaging  is 
used  to  view,  or  image,  an  object  whose  characteristics  of  interest  are  spatial  phase  variations.  Schlieren 
imaging  is  accomplished  here  by  using  a  spatial  filter  in  the  Fourier  plane  of  the  AO  cell  of  interest.  An 
example  of  a  schlieren  image  is  shown  in  Figure  4-5.  This  figure  shows  the  schlieren  image  of  the  Newport 
Electro-Optical  Systems  (NEOS)  multichannel  AO  cell  with  an  RF  signal  applied  to  four  adjacent  channels. 
Through  this  imaging  technique,  AO  cell  acoustic  apertures  can  be  analyzed  and  characterized  to  assess 
such  acoustic  column  properties  as  spreading,  crosstalk,  and  time  aperture. 

Coherent  operation  is  achieved  with  the  in-line  architecture  by  adding  a  reference  oscillator  to  the 
AO  cell  input  to  create  a  reference  beam  to  spatially  and  temporally  interfere  with  the  information  beam.  The 
interference  between  the  reference  beam  and  the  information  beam  at  a  relative  angle  will  create  the  spatial 
fringe  pattern  described  above.  In  addition,  due  to  the  different  Doppler  shifts  induced  by  the  AO  cell,  the 
fringe  pattern  will  oscillate  at  the  difference  frequency  between  the  reference  frequency  and  the  information 
frequency.  Modulation  of  the  laser  diode  cancels  this  beat  frequency  in  order  to  freeze  the  fringe  pattern 
and  allow  for  time  integration.  The  undiffracted  beam  can  be  used  to  interfere  with  the  information  beam, 
but  the  resulting  spatial  carrier  frequency  will  be  prohibitively  high,  requiring  a  very  long  detector  array  to 
adequately  sample  the  carrier.  Using  the  reference  oscillator  allows  the  user  to  set  the  output  carrier  at  the 
desired  frequency.  The  condition  on  this  carrier  is  that  the  reference  oscillator  be  at  least  1 .5  times  the 
bandwidth  away  from  the  information  center  frequency.  Thus,  for  an  80-MHz  center  and  a  10-MHz 
bandwidth,  the  reference  oscillator  must  be  at  a  frequency  of  65  MHz  or  less  (or  greater  than  95  MHz).  For 
this  numerical  example,  with  a  65-MHz  reference  oscillator  resulting  in  a  15-MHz  oscillator  offset,  the  laser 
diode  must  be  driven  with  the  10-MHz  bandwidth  signal  on  a  15-MHz  carrier. 

Mathematically,  these  operations  are  summarized  as  follows.  The  laser  diode  intensity 
modulation  is  given  by. 

I,(t)  =  A^  (1 -t-m,Re{Ag(t)exp(j(t)g(t))exp[j27tfQt)}),  (4-5) 

where  m^  is  the  modulation  depth.  Each  channel  of  the  AO  ceil  is  amplitude  modulated  to  provide  the  spatial 
and  temporal  signal; 

A2(t-t)  =  7^{1 m2A^(t-x)exp[-j<|)|^(t-T)]exp[-j2nfQ(t-T)] } 

•  exp[-j2jtL(t-x)] 

^  '  c'  (4.6) 

where  m2  is  the  ratio  of  signal  amplitude  to  reference  oscillator  amplitude.  Taking  the  squared  magnitude 
of  this  amplitude  and  multiplying  by  the  laser  intensity  modulation  given  above  provides  the  detected 
intensity  for  the  i*^  correlation. 


NEOS  MULTICHANNEL  AO  CELL 
N43080-8-BO 
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Figure  4-5.  Schlieren  Image  of  NEOS  AO  Cell 


ljj(t:T)  =  A^A2(1  +m| 


A,g(t-T)exp[-j(t>^(t-T)]|^) 

2I  2 

+  A^A2m,  (1  +  m||  Af^(t-x)exp(-j4>j^(t-T)]  ) 

•  Re  { Ag  (t)  exp  (t)  ]  exp  [  j2jtfQt] ) 

+  A^  A22m2Re  { A^  (t  -  x)  exp  (t  -  x)  ]  exp  [-j2jtfQ  (t  -  x)  ] } 

(4-7) 

+  A^  A2m^  mgRe  { A^  (t)  exp  (t)  ] 

•  A,g  (t  -  X)  exp  [-j4>^  (t  -  X)  ]  exp  (j2nfQX]  } 

+  A^  Agm^  rrigRe  { A  ^  (t)  exp  (j4>g  (t)  1  A^j  (t  -  x) 

•  exp  (t  -  X)  1  exp  lj2ji  (2fQt  -  IqX)  ]  } . 

The  second,  third,  and  fifth  terms  shown  have  a  temporal  modulation  at  frequency  fo  and  will  integrate  to 
zero  over  typical  detector  integration  times.  The  first  term  is  low  frequency  in  both  space  and  time  relative 
to  the  desired  term,  which  simply  results  in  an  intensity  bias.  The  fourth  term  in  this  equation  is  the  desired 
term  relative  to  the  desired  term,  which  when  integrated  over  time  yields  the  carrier-modulated  correlation 
given  previously  in  Equations  4-3  and  4-4.  After  detector  readout,  but  before  digitization,  it  is  desirable  to 
employ  a  band-pass  filter  (BPF)  to  eliminate  both  the  tow-frequency  bias  and  the  high-frequency  clock 
noise. 

The  detector  sampling  requirement  is  calculated  according  to  standard  Nyquist  sampling  theory. 
First  consider  the  imaging  of  the  S-ps  correlation  window  across  the  linear  array.  The  highest  temporal 
frequency  (encoded  spatially)  for  the  numerical  example  above  is  20  MHz  (15-MHz  center  frequency  with 
a  10-MHz  bandwidth).  One  cycle  at  this  frequency  corresponds  to  50  ns,  which  implies  a  required  pixel 
spacing  of  25  nsec.  Using  this  spacing  over  a  5-ps  window  requires  a  total  of  200  pixels.  Thus,  the  Nyquist 
sampling  requirement  will  be  satisfied  with  a  standard  512-pixel  array.  If  the  carrier  frequency  is  increased 
to  greater  than  15  MHz  (by  using  a  reference  oscillator  of  less  than  65  MHz),  more  than  200  pixels  will  be 
required.  This  shows  why  the  use  of  the  undiffracted  beam  as  a  reference  oscillator  is  impractical  (the 
reference  oscillator  offset  is  80  MHz  requiring  more  than  800  pixels). 

Tests  were  performed  to  determine  the  dynamic  range  of  the  laser  diode  to  be  used  in  the  in-line 
correlator  initial  breadboard.  A  DO  Industries  GALA  laser  system  (model  067-02-04,  670  nm,  5  mW 
maximum  power)  was  chosen  as  the  laser  diode.  A  high-speed  photodiode  (Antel  Corporation)  was  used 
to  detect  the  temporal  modulation  of  the  laser  diode.  The  output  of  the  photodiode  was  observed  on  a 
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Tektronix  11403  digitizing  oscilloscope.  The  dynamic  range  was  calculated  in  the  following  manner.  A 
sinusoid  was  input  to  the  laser  diode.  The  signal  power  was  set  to  correspond  to  the  maximum  output  of 
the  laser  diode.  After  photodiode  detection,  the  input  signal  to  the  laser  diode  was  attenuated  and  the 
corresponding  results  observed  on  the  oscilloscope.  The  amount  of  attenuation  that  still  resulted  in  the 
detection  of  the  modulated  waveform  was  used  as  an  indication  of  the  dynamic  range.  For  the  test,  a 
70  MHz  sinusoid  was  input  into  the  laser  diode.  After  phone  conversations  with  DO  Industries,  it  was 
determined  that  an  upper  limit  on  input  signal  voltage  is  15  Vpp.  It  was  decided  to  use  a  10  Vpp  input  signal 
to  the  laser  diode,  which  provides  a  90%  modulation  of  the  optical  beam.  Figure  4-6  (a)  illustrates  the 
oscilloscope  trace  for  the  1 0  Vpp  input  signal  after  photodiode  detection.  After  46  dB  of  electrical 
attenuation,  a  70  MHz  sinusoid  was  still  visible  on  the  oscilloscope,  as  shown  in  Figure  4-6  (b).  As  noted 
in  the  figure,  the  detected  attenuated  signal  is  very  noisy,  possibly  due  to  detector  noise.  As  a  result  of  the 
tests,  it  was  determined  that  the  dynamic  range  of  the  laser  diode  was  on  the  order  of  46  dB.  This  will  be 
adequate  for  the  MADOP  demonstration  program. 

A  single-channel,  in-line  correlator  was  designed,  constructed,  and  tested  at  the  Photonics 
Center  with  laboratory-generated  signals.  As  mentioned  previously,  an  extension  of  this  design  to  a 
multichannel,  in-line  correlator  would  require  a  multichannel  AO  cell  and  additional  beam  replication  optics. 
A  conceptual  system  layout  for  multichannel  operation  was  shown  in  Figure  4-4.  A  more  detailed  layout  of 
the  single-channel,  in-line  correlator  design  is  illustrated  in  Figure  4-7  (a).  The  combination  of  the  45° 
prism  line  generator  and  piano-cylindrical  lens  CY1  results  in  a  collimated  laser  line  used  to  illuminate  the 
AO  cell.  Spherical  lenses  Si  and  S2  are  used  to  image  the  AO  cell  onto  a  linear  CCD  array.  Note  that  the 
CCD  array  and  the  AO  cell  used  in  this  setup  could  be  replaced  by  the  multichannel  AO  cell  and  the  two- 
channel  line-scan  camera  system  described  in  Subsection  4.1.  The  AO  material  for  this  IntraAction  AOD-70 
device  was  flint  glass  (n=1 .68  @  670  nm,  thickness  ~  60  mm).  The  optical  path-length  due  to  the  index  of 
refraction  and  thickness  of  the  AO  cell  was  taken  into  consideration  when  setting  up  the  imaging  system 
consisting  of  Si  and  S2.  During  setup,  it  was  assumed  that  the  acoustic  column  object  was  located  at  the 
center  of  the  crystal. 

Figure  4-7  (b)  illustrates  the  design  of  the  electronics  used  to  generate  the  input  signals  to  the 
laser  diode  and  to  the  AO  cell.  Note  that  the  dc  bias  term  for  the  laser  diode  input  is  generated  internal  to 
the  GALA  laser  system.  Also,  the  electronics  design  shown  in  the  figure  is  configured  for  realizing  an 
autocorrelation.  Thus,  s{t)  represents  the  temporal  signal  to  be  autocorrelated.  As  shown  in  the  figure,  the 
offset  frequency,  fo.  was  chosen  to  be  10  MHz.  The  requirement  on  the  tone  separation  depends  on  the 
bandwidth  of  information  being  input  to  the  correlator  and  is  given  by: 

fo>l.5B  (4-8) 
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-19.8  ns  lOns/div  80.2  ns 

(b)  s  (t)  =  70  MHz  Sinusoid.  After  46  dB  Attenuation 
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Figure  4-6.  Laser  Diode  Dynamic  Range  Measurements  (Concluded) 
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(b)  Electronic  Signal  Generator  Design 


Figure  4-7.  In-Line  Correlator  Breadboard  Design 


4-14 


where  Iq  is  the  frequency  separation  between  the  IF  of  the  information  and  the  reference  tone  frequency, 
and  B  is  the  largest  bandwidth  of  information  input  into  either  the  laser  diode  or  AO  cell.  The  upper  limit  on 
fo  depends  on  the  bandwidths  of  the  AO  cell  and  the  laser  diode,  and  the  sampling  of  the  CCD  array. 

Figure  4-8  (a)  illustrates  the  theoretical  prediction  of  the  autocorrelation  of  a  pulse  train  whose 
parameters  are  listed  in  the  figure.  Given  a  pulsewidth  of  t  and  a  pulse-to-pulse  separation  of  T  (T=1/prf), 
the  expected  triangular  correlation  theoretically  has  a  width  of  2t  and  a  correlation  peak-to-peak  separation 
of  T.  Thus,  for  a  pulse  train  that  has  a  pulsewidth  of  t  =  0.5  ps  and  T  =  2.5  ps,  the  expected  autocorrelation 
width  would  be  1  ps  with  correlation  separation  of  2.5  ps.  Since  5.0  ps  of  the  AO  cell  aperture  was  imaged 
across  the  CCD,  it  was  expected  that  there  would  be  two  correlation  peaks  across  the  CCD  array. 
Figure  4-8  (b)  shows  the  autocorrelation  of  s(t)  at  the  output  of  the  CCD  array.  Note  that  another  detector 
with  a  1 2-bit  A/D  conversion  system  was  employed  during  this  data  collection.  The  spatial  carrier  of  the 
correlation  corresponds  to  a  reference  tone  separation  ot  10  MHz,  and  the  integration  time  for  the  CCD 
array  was  nominally  set  at  3  ms.  As  seen  in  Figure  4-8  (b),  the  experimental  results  were  very  similar  to 
the  theoretical  results  described  earlier.  The  additional  background  ripple  is  a  result  of  an  additional  bias 
term  imparted  on  the  signal,  s(t). 

Figure  4-9  (a)  shows  s(t)  for  a  pulse  train  for  which  t  =  0.25  ps  and  T  =  1.25  ps.  Note  the 
additional  bias  term  previously  mentioned  results  in  the  background  ripple.  The  correlation  should  be  similar 
in  shape  to  that  in  Figure  4-8  (b).  However,  it  is  expected  that  the  correlation  width  and  separation  would 
be  half  of  that  given  in  the  previous  example  and  that  the  number  of  correlation  peaks  would  double  over 
the  5.0  ps  imaged  time  aperture.  Figure  4-9  (b)  shows  the  autocc  relation  of  s(t)  at  the  output  of  the  CCD 
array  and  the  agreement  with  expected  results.  As  in  the  previous  example,  the  spatial  carrier  of  the 
correlation  corresponds  to  a  reference  tone  separation  of  10  MHz  and  the  integration  time  for  the  CCD  array 
was  set  at  3  ms. 

As  mentioned  earlier,  the  spatial  carrier  for  the  in-line  architecture  can  be  modified  by  changing 
the  separation  of  the  electronic  reference  tone  with  respect  to  the  information  at  the  IF.  Increasing  the 
separation  between  the  reference  tone  and  the  IF  increases  the  spatial  frequency  of  the  correlation  carrier. 
Figure  4-9  (c)  shows  the  autocorrelation  of  s(t)  for  a  pulse  train  for  which  t  =  0.25  ps,  T  =  1 .25  ps,  and  fo  = 
15  MHz.  Comparing  with  Figure  4-9  (b).  it  is  obvious  that  the  spatial  carrier  frequency  has  increased  for 
increasing  tone  separation.  There  are  restrictions/limits  that  must  be  taken  into  consideration  when 
selecting  the  tone  separation.  As  mentioned  earlier,  the  minimum  separation  is  set  by  the  bandwidths  to  be 
correlated  (see  equation  4-8).  The  upper  limit  on  the  tone  separation  is  a  function  of  the  number  of  pixels 
in  the  CCD  array  and  the  total  delay  window.  There  must  be  an  adequate  number  of  pixels  to  sample  each 
fringe.  Other  signals  were  tested  with  the  in-line  correlator,  including  square  waves  and  sinusoids.  These 
results  are  documented  in  the  Dynetics  on-site  lab  notes. 
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Figure  4-8.  Autocorrelation  of  »(t)  (In-Line):  eft)  =  Pulse,  T=  0.5  ns, 
T^2.5^f0=  10  MHz 
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Figure  4-9.  In-Line  Autocorrelation  Results  (Concluded) 
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Preliminary  testing  was  performed  to  determine  an  estimate  of  the  in-line  correlator  dynamic 
range.  The  method  used  was  similar  to  that  for  testing  the  dynamic  range  of  the  laser  diode.  A  pulse  signal, 
s(t),  with  pulsewidth  t  =  1 .0  ps,  pulse  repetition  interval  T  =  5.0  ps,  and  peak-to-peak  amplitude  A  =  1  Vpp, 
was  input  to  the  in-line  correlator  system.  Figure  4‘10  (a)  illustrates  the  autocorrelation  of  the  pulse 
waveform  with  the  above  parameters.  The  waveform  was  attenuated  prior  to  insertion  into  the  laser  diode 
(the  power  into  the  AO  cell  was  kept  constan-).  This  would  be  similar  to  the  reduction  of  the  error  signal  in 
the  MADOP  demonstration.  Figure  4-10  (b)  and  4-10  (c)  shows  the  autocorrelation  of  s(t)  after  25  and  33 
dB  of  signal  attenuation.  Even  after  33  dB  of  attenuation,  the  correlation  of  the  signal  still  can  be  observed. 
During  this  set  of  tests,  however,  the  background  nonuniformities  were  significant  on  the  CCD  array.  The 
background  nonuniformities  can  be  seen  by  inputting  a  single  tone  to  the  AO  cell  and  switching  off  the 
modulation  to  the  laser  diode.  The  resultant  CCD  output  is  illustrated  in  Figure  4-11.  Possible  sources  for 
this  non-uniformity  may  include  acoustic  attenuation  in  the  AO  cell  and  illumination  beam  nonuniformities. 
Reduction  of  this  noise  can  be  accomplished  through  post-processing  techniques,  as  will  be  illustrated  in 
the  next  section  for  the  two-path  correlator  architecture.  At  the  time  that  these  data  were  recorded,  no  post¬ 
processing  software  had  been  developed. 

The  testing  performed  on  the  in-line  correlator  has  given  much  insight  to  the  architecture  setup, 
operation,  capabilities  and  limitations.  Further  analysis/testing  should  be  performed  since  operation  and 
performance  of  the  new  components,  such  as  the  multichannel  AO  cell  and  CCD  array,  may  differ  from  that 
of  the  current  components  in  the  preliminary  setup.  Tests  that  should  be  performed  include  further  dynamic 
range  characterization,  linearity,  delay  window  uniformity,  and  the  effects  of  the  post-detection  processing 
algorithms  on  these  measurements.  Due  to  time  and  hardware  constraints,  full  characterization  of  the  in¬ 
line  architecture  did  not  occur. 

4.2.2  Design  Alternative  II:  Two-Path  Multichannel  Correlator  Architecture 

An  alternative  to  the  above  in-line  design  is  achieved  through  the  application  of  two  AO  cells  in 
independent  paths  of  a  Mach-Zehnder  interferometer,  as  shown  in  Figure  4-12.  A  continuous  wave  (CW) 
laser  source  is  split  into  two  separate  paths  through  the  use  of  a  beam  splitter.  Both  paths  of  the 
interferometer  contain  AO  cells  through  which  the  auxiliary  antenna  signals  and  the  error  signal  will  be  input 
to  the  correlator.  The  error  signal  is  input  to  the  single-channel  AO  cell  and  the  auxiliary  antenna  inputs  are 
connected  to  separate  channels  of  a  multichannel  device.  The  error  channel  can  be  replicated  in  the  same 
manner  as  the  laser  diode  signal  replication  described  previously,  or  through  the  use  of  another 
multichannel  AO  cell  in  the  error-signal  path.  In  the  latter  case,  additional  RF  components  would  be  needed 
to  electrically  replicate  the  RF  signal  and  input  the  signals  into  separate  channels  of  the  second 
multichannel  device.  Advantages  of  using  another  multichannel  AO  cell  in  the  error  channel  path  will  be 
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(a)  Autocorrelation  of  s(t)  =  Pulse,  t  =  1  |is,  T  =  5  ^s,  A  =  1  Vpp 


(b)  Autocorrelation  after  25  dB  Attenuation  of  s(t) 
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Figure  4-10.  In-Line  Correlator  Dynamic  Range  Measurements 
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GRAYSCALE  VALUE 


(c)  Autocorrelation  after  33  dB  Attenuation  of  s(t) 
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Figure  4-10.  In-Line  Correlator  Dynamic  Range  Measurements  (Concluded) 


4-22 


GRAYSCALE  VALUE 


PIXEL  NUMBER 

TR-92-SI0 1-0436 


Figure  4-11.  Illumination  of  CCD  Array  •  No  Modulation  Input  to  Laser 
Diode,  Tone  into  AO  Cell 
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described  later.  The  separate  paths  are  recombined  through  another  beamsplitter  and  optical  components 
are  used  to  schlieren  image  the  AO  cell  diffracted  beams  onto  linear  CCD  arrays. 

In  this  architecture,  the  spatial  carrier  frequency  required  to  operate  coherently  is  set  by  tilting  the 
beam  combiner  cube  to  offset  the  images  on  the  detector  array  in  angle.  Because  the  recombined  beams 
from  a  given  channel  have  been  equally  Doppler-shifted  by  the  AO  interaction  (the  same  diffraction  order 
must  be  selected  for  each  path),  the  correlation  fringe  pattern  is  stationary,  thus  allowing  CW  illumination  of 
the  AO  cells. 

Mathematically,  this  architecture  can  be  viewed  as  follows.  The  amplitude  of  the  single-channel 
AO  cell  can  be  expressed  as; 

A^  (t  -  X)  =  Ag  (t  -  X)  exp  [-j<t)g  (t  -  X)  ]  exp  [-j2nfQ  (t  -  x)  ] ,  (4-9) 

and  for  a  channel  of  the  multichannel  AO  ceil  as; 

A2(t  +  x)  =  A,^j(t  +  x)exp(-j<>^(t  +  x)]exp[-j2jrfQ(t-i-x)] .  (4-10) 

It  can  be  seen  from  these  equations  that  the  propagation  directions  are  counterpropagating.  The  detector 
sees  the  squared  magnitude  of  the  sum  of  these  two  amplitudes  resulting  in  the  detected  signal; 

'det  =  I  \  ®^P  H*!*®  ~  1  H2«fQ  (t  +  X)  ] 

+  A,^(t  +  x)expt-]())j^(t  +  x)]expI-j27rfo(t  +  x))|2 

=  |A^(t-x)|2  +  |A,^(t  +  x)|2 

+  2Re  { A  ^  (t  -  X)  exp  [- (t  - x)  ]  A^  (t  +  x) 


•exp[j<|>^(t  +  x)]exp[j4nfQX]  }. 

where  ldet(f  !'<)  'S  the  instantaneous  intensity  at  the  detector.  Letting  t'  =t-x,  this  equation  as  can  be  written  as; 
ldgt(t';^)  =  low  frequency  biases 

+  2Re  J  Ag(n  expH<{)^(t')]  A,^(r  +  2x) 

0 

•  exp  [j<t>^  (f  +  2x)  1  exp  [j4nfQx]  dt'j 

(4-12) 
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It  can  be  seen  from  this  coordinate  transformation  that  the  effective  correlation  window  is  doubled  for  the 
counterpropagating  signals  relative  to  the  noncounterpropagating  situation  characteristic  of  the  in-line 
alternative.  Here  again,  only  one  term  is  at  a  spatial  carrier  and  temporally  at  dc.  Integration  of  this  term 
results  in  the  desired  carrier-modulated  correlation  given  above.  The  spatial  carrier  frequency  shown  above 
is  the  result  of  the  diffraction  angle  from  each  cell  (a  function  of  the  center  frequency),  and  as  previously 
mentioned,  can  be  modified  to  the  desired  spatial  carrier  through  adjustment  of  the  beam  combiner  cube. 

In  this  design  of  the  two-path  architecture,  the  error  signal  is  input  to  a  single-channel  AO  cell  and 
the  auxiliary  returns  are  input  to  the  multichannel  AO  cell.  The  diffracted  beam  from  the  single-channel  AO 
cell  is  replicated  through  the  use  of  beam  splitting  optics  for  interference  with  the  multi-channel  diffracted 
beams  at  the  final  beam  combiner.  As  a  result,  the  straight-through,  single-channel  diffracted  beam  would 
interfere  with  one  of  the  multichannel  diffracted  beams,  and  the  replicated  single-channel  beam  would 
interfere  with  the  other. 

A  multichannel,  two-path,  time-integrating  correlator  was  designed,  implemented  in  hardware, 
and  tested  with  laboratory-generated  signals.  The  purpose  of  this  stage  of  design  was  to  determine  the 
utility  of  the  beam  splitter/mirror  assembly  proposed  as  the  beam-replication  optics.  Figure  4-13  shows  the 
lavout  of  the  two-path  correlator  with  the  incorporation  of  the  two  1  x  512  pixel  EG&G  LC1901FKN-011 
modular  line  scan  cameras  and  the  display  of  the  correlations  through  the  PC  and  the  MAC  computer 
systems.  Software  developed  by  1 LT  Andrews  allowed  for  real-time  display  of  the  camera  outputs.  With 
the  linear  CCD  array  and  display  software  previously  utilized,  real-time  display  of  the  CCD  array  output  was 
not  possible. 

The  single-channel  AO  cell  used  in  this  architecture  is  identical  to  the  AO  cell  used  in  the  in-line 
hardware  implementation.  The  multichannel  AO  cell  used  for  the  two-channel  test  setup  is  a  NEOS  device. 
The  index  of  refraction  and  thickness  of  the  AO  cells  and  the  beam  splitter  (part  of  the  beam  replicating 
optics  in  the  error  signal  path)  are  given  in  Table  4-1.  This  information  was  used  to  calculate  the  path 
lengths  due  to  the  above  components,  interferometer  path  length  differences,  and  imaging  compensation 
due  to  these  components.  Imaging  compensation  refers  to  the  ^act  that  there  is  a  common  image/ 
interference  plane  and  that  AO  cell  locations  from  the  imaging  lenses  are  adjusted  accordingly. 

Unlike  the  in-line  architecture,  where  the  electronic  reference  determines  the  spatial  carrier 
frequency,  the  spatial  carrier  in  the  two-path  architecture  is  set  by  the  angles  at  which  the  beams  in  the  two 
paths  interfere.  This  angle  is  controlled  by  the  tilt  and  rotation  of  the  second  beam  splitter/combiner,  (BS2). 
Referring  to  Figure  4-14,  rotation  of  BS2  in  the  y  direction  results  in  changing  the  tilt  angle  of  fringes.  Also, 
rotation  of  BS2  in  the  x  direction  results  in  changing  the  frequency  spacing  of  the  spatial  fringes.  Note  that 
in  both  cases,  the  correlation  envelope  width  stays  the  same  and  only  the  tilt  and  rotation  of  the  fringes 
within  the  correlation  envelope  changes.  Wider  separation  of  fringes,  accomplished  by  ±  x  direction 
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Figure  4-13.  Two-Path  Time-Integrating  Correlator 
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movement,  leads  to  a  "zero  fringe"  condition  during  which  the  correlation  peak  appears  to  be  either  all  light 
or  all  dark.  This  is  illustrated  in  Figure  4-14  (c).  Continuing  to  rtxjve  in  the  x  direction  leads  to  increasing 
fringe  frequency  (smaller  fringe  separation).  From  past  experience,  it  is  highly  desirable  to  have  vertical 
fringes  if  the  CCD  pixels  and  acoustic  apertures  are  oriented  along  the  horizontal  axis.  If  a  cylindrical  lens 
is  used  to  collapse  the  information  onto  a  linear  CCD  array,  vertical  fringes  will  aid  in  increasing  fringe 
contrast  by  minimizing  collapsed  fringe  overlap  on  the  CCD  array. 


Table  4-1.  Optical  Path  Length  Parameters 


Component 

Index 
at  633  nm 

Thk':kness 

(mm) 

Single-Channel 
AO  Cell 

1.68  (flint  glass) 

6CI.0 

Multichannel 

AO  Cell 

2.26  (Te02) 

20,0 

Beam  splitter 

1.515  (glass) 

50.8 

Multichannel  correlation  was  achieved  using  a  5  MHz  noise  source  (centered  at  the  system  IF) 
as  the  system  input.  A  1  MHz  square-wave  input  signal  was  DSB-SC  modulated  at  80  MHz  and  input  to 
the  AO  cells  as  a  second  test  waveform.  Figures  4-15  and  4-16  show  the  multichannel  correlation  results 
as  displayed  on  the  PC  monitor  and  from  the  2-b  CCD  array. 

In  Figures  4-15  (b)  and  4-16  (b),  the  correlation  pattern  toward  the  bottom  of  both  figures  is  that 
from  the  replicated  signal  path.  As  can  be  seen  in  the  figures,  there  appear  to  be  some  additional  distortions 
in  the  replicated  path.  This  can  be  best  visualized  in  Figure  4-16  (b)  by  noting  the  "S”  pattern  of  the  fringes 
for  the  replicated-path  correlation  pattern.  This  was  the  best  fringe  pattern  obtainable  through  the  use  of 
the  beam  splitter/mirror  assembly  adjustments.  One  cause  of  the  distortions  may  be  the  surface  quality  of 
the  flat  mirror.  With  this  in  mind,  it  was  decided  that  in  the  MADOP  hardware  implementation,  it  would  be 
easier,  from  an  alignment  standpoint,  to  use  a  multichannel  AO  cell  in  the  error  path  of  a  two-path 
architecture.  The  insertion  of  the  second  multichannel  AO  cell  into  the  error  path  of  the  correlator  would 
allow  easy  replication  of  the  error  channel  by  signal  splitting  and  inputting  the  error  signal  into  two  separate 
channels  of  the  multichannef  device.  Additional  advantages  of  using  a  multichannel  AO  cell  in  the  error  path 
of  the  correlator  include  the  elimination  of  the  beam  splitter/mirror  assembly,  the  reduction  of  optical  path 
differences  between  correlator  paths  and  between  correlator  channels,  and  co-planar  diffracted  beams  from 
the  AO  cell  for  each  channel  utilized  (reducing  alignment  requirements  of  channels  in  the  error  path).  It  is 
important  that  the  fransducer/channel  spacing  is  the  same  for  both  the  multichannel  cells,  since  it  is  very 
difficult  to  adjust  the  individual  channel  overlap  in  the  correlation/acoustic  image  plane.  Thus,  it  is  desirable 
to  procure  an  additional  matched  (in  terms  of  device  specifications  and  operational  parameters) 
multichannel  AO  cell  to  be  used  as  the  error  signal  AO  cell . 
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Figure  4-14.  Spatial  Fringe  Frequency/Orientation  as  Function 
of  Acoustic  Beam  Image  Overlap  Positions 


(a)  Linear  CCD  Array  Output 
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Figure  4-15.  Multichannel  Correlation  of  5  MHz  Noise  Source  Input 
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(b)  2-D  CCD  Output 

Figure  4-15.  Multichannel  Correlation  of  5  MHz  Noise  Source  Input  (Concluded) 


(a)  Linear  CCD  Array  Output 
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Figure  4-16.  Multichannel  Correlation  of  1  MHz  Square  Wave  Input 
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(b)  2-D  CCD  Output 

Figure  4-16.  Multichannel  Correlation  ofl  MHz  Square  Wave  Input  (Concluded) 


A  redesign  of  the  two-path  correlator  optical  layout  was  performed  due  to  the  insertion  of  the 
Brimrose  multichannel  AO  cell  (in  the  auxiliary  antenna  path)  and  the  use  of  the  EG&G  CCD  arrays. 
Figure  4-77shows  a  top-view  optical  layout  of  the  interferometric  paths  of  the  two-path  correlator.  L1  and 
a  second  lens  L2  (not  shown),  placed  after  BS2,  are  used  to  image  the  appropriate  AO  cell  apertures  back 
onto  the  CCD  arrays.  The  magnification  of  the  imaging  system  is  given  by  the  ratio  of  the  focal  lengths  of 
L1  and  L2.  Figure  4-7 7 illustrates  how  path  differences  due  to  different  AO  cell  materials  and  thicknesses 
were  taken  into  consideration  during  component  layout  and  correlator  alignment.  In  the  figure,  the  two 
optical  paths  of  the  Mach-Zehnder  interferometer  are  laid  out  without  path  folding.  The  top  path  is  the  one 
with  the  interaction  AO  cell  (A01)  and  the  bottom  path  is  the  Brimrose  AO  cell  (A02)  path.  BS1  and  BS2 
refer  to  the  first  and  second  beam  splitters,  respectively:  t^  and  t2  refer  to  the  thicknesses  of  the  AO  cells: 
ni  and  n2  are  the  indices  of  refraction  of  the  AO  cells:  SF  refers  to  the  spatial  filters:  A  is  a  set  separation 
between  SF  and  BS2:  and  and  A2  refer  to  the  distance  between  BSl  and  the  respective  AO  cell.  As 
indicated,  L2  is  located  after  the  second  beam  splitter  (BS2).  For  the  optical  component  layout  calculations, 
the  center  of  the  acoustic  column  location  was  assumed  to  be  located  at  t/2.  This  gives  rise  to  the  nj  \J2 
location  of  the  object  plane  from  the  front  face  of  the  AO  cell.  The  total  path  length  in  each  path  is  given  by 
the  following  equation: 


L(1)  =  2f^  +A  +  A^ 


’4) 


A01 


L(2)  =  2f^  +  A-t- A^-i- 


(4-13) 


where  L(i)  is  the  total  path  length  in  path  i  from  BSl  to  8S2,  and  f  1  is  the  focal  length  of  LI . 

The  requirement  for  optical  system  magnification  was  calculated  in  the  following  manner.  The 
Brimrose  AO  cell  is  Te02  [L]  and  has  an  acoustic  velocity  of  4.2  X  1 0^  m/s.  Thus,  a  2.5  ps  time  aperture  in 
the  AO  cell  corresponds  to  approximately  10.5  mm.  In  the  IntraAction  AO  cell,  this  aperture  dimension  was 
calculated  to  be  10.05  mm.  Note  that  with  the  two-path,  Mach-Zehnder,  time-integrating  correlator 
architecture,  the  correlation  time  window  is  twice  that  of  the  imaged  AO  cell  apertures  if  equal  acoustic 
apertures  are  imaged  (see  Equation  4-12).  Given  a  6.656-mm  detector  array  length  for  the  EG&G  CCD 
arrays,  a  magnification  of  M  =  0.634  was  calculated  to  image  10.5  mm  of  each  AO  cell.  Using  a  two-lens 
imaging  system  and  the  available  lenses,  a  magnification  of  M  =  0.625  was  achieved  using  f,  =  200  mm 
and  f2  =  125  mm.  With  these  focal  lengths,  the  path  length  for  PdCh  leg  of  the  interferometer  was  selected 
to  be  L(1)  =  L(2)  =  488  mm. 

Figure  4-78  illustrates  the  design  of  the  electronics  used  to  generate  the  input  signals  to  the  AO 
cells.  The  electronics  were  designed  to  generate  DSB-SC  signals  at  the  system  IF  (80  MHz).  Test  signals 
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Figure  4-17.  Interferometic  Path  Optical  Layout 


were  generated  at  base  band  and  mixed  with  a  carrier  tone  of  80  MHz.  In  the  figure,  I,  L,  and  R  refer  to  the 
input/output  nomenclature  on  the  Mini  Circuits  Labs  (MCL)  mixer.  It  should  be  noted  here  for  reference  that 
when  using  MCL  mixer  P/N  ZP-10514,  DSB-SC  could  not  be  achieved.  Switching  to  an  MCL  mixer  P/N 
15542  (as  noted  in  the  figure)  enabled  DSB-SC  signal  modulations  to  be  accomplished.  In  addition  to  using 
band  limited  (at  base  band)  noise  signal  inputs,  the  LeCroy  Arbitrary  Function  Generator  (AFG)  was  used 
to  generate  some  base  band  test  signals  which  were  DSB-SC  modulated  prior  to  system  input. 
Figure  4-19  (a)  and  (b)  illustrates  the  autocorrelation  of  a  500  kHz  sinusoid  as  viewed  by  the  EG&G  linear 
CCD  array  and  by  the  COHU  2-D  CCD  arrays,  respectively. 
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Figure  4-18.  Electronics  Design  for  Signal  Inputs  to  AO  Cells 


Enhancements  to  the  PC  display  software,  implemented  by  1LT  Andrews,  increased  the 
capabilities  of  the  real-time  EG&G  array  displays.  For  example,  menu-driven  options  exist  for  the 
subtraction  of  array  backgrounds.  This  option,  along  with  the  addition  of  a  positive  bias,  was  utilized  for  the 
subtraction  of  beam  nonuniformities  due  to  the  illumination  source  and  due  to  the  AO  cell  diffraction  pattern. 
Implementation  of  this  was  performed  in  the  following  manner.  The  signals  were  input  to  the  AO  cells  in  the 
regular  manner  for  detection  of  autocorrelation.  This  included  adjusting  the  power  level  into  the  AO  cells 
so,  on  the  average,  the  diffracted  intensities  from  the  two  cells  were  equal.  With  one  of  the  optical  paths 
blocked,  the  intensity  pattern  incident  on  the  CCD  array  showed  the  nonuniformities  described  above.  After 
subtracting  the  background  pattern  (from  a  preset  positive  bias),  the  above  step  was  repeated  for  the  other 
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(b)  COHU  CCD  Array  Output 

Figure  4-19.  Correlation  Results;  s  (t)  =  500  kHz  Sinusoid  (Concluded) 
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optical  path.  This  resulted  in  a  uniform  bias  about  which  the  correlation  signal  was  positioned. 
Figure  4-20  (a)  illustrates  the  correlation  output  for  wide  band  noise  input  to  the  AO  cells.  The  overlapping 
3  dB  bandwidth  of  the  signal,  set  by  the  AO  cell  pass  bands,  was  25  MHz.  As  shown  in  the  figure,  there  is 
a  nonuniform  background/bias.  Figure  4-20  (b)  illustrates  the  same  correlation  output  after  bias  addition 
and  subtraction  of  the  background  in  the  manner  described  above.  Note  the  improvement  in  the 
background  signal  and  the  increased  correlation  peak  modulation  depth.  There  are  other  methods  for 
performing  this  sort  of  "post-processing”  technique,  such  as  bi-phase  substraction  and  bandpass  filtering, 
but  for  the  current  stage  of  correlation  testing,  the  current  near-real-time  technique  is  adequate. 

A  preliminary  test  of  the  dynamic  range  of  the  two-path  correlator  was  performed.  A  pulse  train 
(t  =  0.1  ps,  T  =  5.0  ps),  shown  in  Figure  4-21  (a),  was  DSB-SC  modulated  at  80  MHz  and  input  into  both 
AO  cells.  Note  that  there  is  no  bias  on  the  pulse  train.  As  mentioned  earlier,  if  the  bias  is  part  of  the  signal 
being  autocorrelated,  the  autocorrelation  of  the  bias  will  also  be  modulated  by  the  carrier.  With  the  bias 
removed  from  the  pulse  signal,  fringes  will  only  be  present  across  the  autocorrelation  of  the  pulse. 
Figure  4-21  (b)  illustrates  the  input  spectrum  to  the  AO  cells  and  Figure  4-21  (c)  shows  the  output 
correlation  pattern  after  performing  the  post-processing  operation  described  above.  The  error  channel  input 
was  attenuated  and  the  correlation  was  viewed  on  the  PC  monitor.  Due  to  the  changing  input  RF  power 
levels  in  the  error  channel,  the  background  subtraction  technique  was  repeated  after  each  attenuation  step. 
The  autocorrelation  of  s(t)  after  30  and  35  dB  attenuation  of  the  error  channel  input  is  shown  in 
Figure  4-21  (d)  and  4-21  (e).  Incorporation  of  additional  or  alternative  post-processing  techniques  may 
increase  this  dynamic  range.  Another  method  for  increasing  the  dynamic  range  would  incorporate  the  use 
of  automatic  gain  control  (AGC)  to  the  error  signal  prior  to  insertion  into  the  AO  cell. 

4.2.3  Time-intearatlna  Correlator  Architectural  Considerations 


It  is  desirable  to  select  and  emphasize  one  of  the  time-integrating  correlator  architectures  with 
respect  to  hardware  fabrication  and  signal  testing.  To  assist  in  the  architectural  evaluations.  Table  4-2  lists 
the  features  associated  with  each  of  the  architectures.  One  of  the  in-line  architecture  drawbacks  that  should 
be  pointed  out  is  that  a  long  time-delay  (on  the  order  of  half  the  desired  time  aperture)  is  required  prior  to 
signal  insertion  into  the  laser  diode.  This  is  necessary  to  center  the  correlation  in  the  delay  window  as 
illustrated  in  Figure  4-22  Consider  a  signal,  s{t),  simultaneously  input  to  a  laser  diode  and  an  AO  cell.  Due 
to  the  nature  of  the  AO  cell,  the  signal  representation  as  it  propagates  through  the  crystal  (neglecting 
propagation  losses)  is  given  by: 


s(t,x) 


X  \ 


(4-14) 
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(a)  Linear  CCD  Array  Output 


! 

. <•■ 


•  1 1  •  !  •  ' 
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(b)  Linear  CCD  Array  Output  After  Background  Subtraction 


Figure  4-20.  Noise  Autocorrelation 


I 
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where  T  is  the  time  aperture  of  the  AO  cell,  x  is  the  dimension  in  the  direction  of  the  acoustic  propagation 
(referenced  to  the  center  of  the  aperture),  and  Vg  is  the  acoustic  velocity  of  the  AO  crystal.  The  term  (T/2  + 
x/Vg]  is  commonly  referred  to  as  the  delay.  When  a  wide  band  signal  is  input  to  the  in-line  correlator  and 
there  is  no  delay  between  the  signal  inputs  to  the  AO  cell  and  the  laser  diode,  the  correlation  occurs  at  the 
zero  delay  point.  As  noted  in  the  figure,  the  zero  delay  point  corresponds  to. 


X  = 


(4-15) 


which  is  at  the  transducer  side  of  the  acoustic  aperture.  To  center  the  correlation  peak  at  the  center  of  the 
aperture,  it  is  necessary  to  time-delay  the  input  to  the  laser  diode  by  an  amount  equal  to  T/2.  For  the  time 
apertures  necessary  for  the  MAOOP,  these  time-delay  differences  can  be  easily  achieved  through  the  use 
of  the  BAW  delay  lines  being  procured.  Another  possible  solution  to  this  issue  involves  using  an  AO  point 
modulator  to  provide  the  necessary  intensity  modulation.  Selection  of  different  delays  can  then  be 
accomplished  through  selecting  different  locations  in  the  rTX)dulator  to  illuminate. 


Table  4-2.  Time-Integrating  Correlator  Architecture  Features 


Two  Path  (Additive) 

In-Line  (Multiplicative) 

Signal  inputs  at  system  IF 

Signal  inputs  to  laser  not  at  system  IF 

Sensitive  to  vibrationAhermal  gradients 

Stability  inherent  in  common-path  architecture 

Manual  spatial  carrier  adjustment 

Electronic  control  of  spatial  carrier 

CW  illumination  of  AO  cells 

Modulated  illumination  source  required  (laser 
diode  or  point  modulator) 

Range  window  is  twice  AO  time  aperture 

Range  window  is  AO  cell  time  aperture 

Cell  bandwidth  -  processing  bandwidth 

Cell  bandwidth  ^1.5  processing  bandwidth 

Signal  generation  of  reference  tones  required 

Vertical  fringes 

Laser  diode  stability  issues 

Large  RF  time  delay  required  in  laser  diode 
path  (half  of  time  aperture).  Needed  to 
achieve  center  correlation 

With  all  the  architectural  considerations  presented  above  and  in  Table  4-2,  it  was  decided  that 
the  emphasis  for  hardware  implementation  and  testing  would  be  the  two-path  correlator  architecture.  It  is 
recommended  that  work  on  the  in-line  correlator  hardware  implementation  proceed  on  an  as-needed  basis. 
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4.3 


AOTDL  FILTER 


Much  of  the  analysis,  design,  and  hardware  implementations  during  the  current  effort  was 
directed  at  the  tapped  delay  line  portion  of  the  system.  The  previous  architecture  (Reference  1)  provided 
effective  tapped  delay  line  filtering,  but  required  DSB  waveform  modulation  together  with  a  strong  tone  at 
the  carrier  frequency.  The  architectures  described  in  this  subsection  are  evolutions  of  this  prior  system  to 
realize  filtering  of  arbitrary  signal  modulations  at  the  system  IF.  Subsection  4.3.1  describes  in  mathematical 
detail  a  Mach-Zehnder  architecture  for  creating  a  tapped  delay-line  that  outputs  the  desired  signal  at  the 
system  80-MHz  IF.  This  architecture  was  found  to  be  deficient  in  generating  true  time-delay  due  to  very 
subtle  properties  of  AO  cells  and  interferometric  architectures.  This  deficiency  will  be  fully  described  in 
Subsection  4.3. 1.7. 

A  noninterferometric  architecture  for  generating  the  desired  tapped  delay-line  properties  was  a 
natural  outgrowth  of  the  increased  insight  into  the  limitations  of  the  Mach-Zehnder  approach.  The 
development  of  a  two-channel  version  of  this  architecture  is  reviewed  in  Subsection  4.3.2.  Included  are 
characterization  results  that  point  to  the  need  for  several  improvements,  most  notably  in  the  stabilization  of 
fhe  system  and  the  optical  lens  design. 
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4.3.1 


The  Mach-Zehnder  architecture  for  implementing  the  tapped  delay  line  is  shown  in  Figure  4-23. 
The  key  to  this  architecture  is  the  addition  of  a  reference  path  that  is  initiated  after  the  AO  spatial  light 
modulator  (AOSLM).  This  path  is  Doppler-shifted  by  the  input  tap  frequency  information  in  an  identical 
fashion  to  the  Doppler-shifted  light  that  illuminates  the  AOTDL.  After  recombination  on  the  photodetector, 
this  effectively  cancels  the  tap-dependent  Doppler  frequency  shift  contribution  of  the  AOSLM,  provides  the 
output  at  the  same  IF  as  the  signal  input  to  the  AOTDL,  and  reduces  the  impact  of  cross-terms  in  this  output. 

Referring  again  to  Figure  4-23.  the  operation  can  be  described  as  follows.  The  AOSLM  is 
illuminated  with  an  unmodulated  "sheet"  beam  that  is  collimated  in  the  direction  of  acoustic  propagation  and 
focused  using  cylindrical  lenses  in  the  orthogonal  direction  (not  shown).  The  output  of  the  AOSLM  is 
deflected  in  angle  proportional  to  the  input  frequencies,  is  Doppler-shifted  by  an  amount  equal  to  these  input 
frequencies,  and  has  an  amplitude  proportional  to  the  amplitude  of  the  input  frequency  components.  Lens 
LI  then  forms  the  spatial  Fourier  transform  of  this  information  at  plane  Pi.  where  an  aperture  allows  only 
the  +1  order  of  diffraction  to  pass.  Lens  L1  has  a  long  focal  length,  f^,  to  provide  a  large  Fourier  transform 
where  the  frc./uericy  coordinate,  t.  is  given  as; 

X  =  XfJ^  =  Xf.  (4-16) 

a 

wrere  X  is  the  optical  wavelength,  f*  is  the  spatial  frequency  of  the  input  acoustic  wave  (f*  =  / /Vg).  /  is  the 
input  frequency  and  Vg  is  the  acoustic  velocity  (4200  m/s  for  longitudinal  Te02).  Lenses  L2  and  L3  then 
serve  as  a  telescope  to  image  this  Fourier  plane,  with  magnification  13/12,  to  the  plane  of  the  AOTDL.  Thus, 
a  single  frequency  at  the  input  to  the  AOSLM  acts  to  tap  the  AOTDL  at  a  given  position: 

x  =  X-^^  (4-17) 

‘2  ''a 

and  with  an  optical  frequency  that  is  Doppler-shifted  by  an  amount  equal  to  /  Lens  L4  then  takes  a  final 
Fourier  transform  of  this  output  to  convert  the  tap  positions  into  tilted  plane  waves  that  are  incident  on  the 
detector  with  an  angle  proportional  to  tap  position.  The  combination  of  lenses  L1  and  L2  (a  telescope) 
images  the  AOSLM  aperture  to  plane  P2,  then  the  combination  of  lenses  L3  and  L4  images  plane  P2  fo  the 
plane  of  the  detector.  Therefore,  fhe  combination  of  AO  cells  and  lenses  serves  fo. 

1 .  Image  the  AOSLM  aperture  onto  the  detector  with  magnification 

M  =  (4-18) 

M'3 
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2.  Transform  the  AOTDL  to  create  plane  waves  incident  on  the  detector  for  each  tap  position  at 
an  angle  proportional  to  the  tap  position,  and 

3.  Transform  the  AOTDL  to  create  images  of  the  AOSLM  at  positions  corresponding  to  each  of 
the  frequencies  input  into  the  AOTDL. 

The  reference  path  is  created  by  a  beam  splitter  at  the  input  to  the  AOTDL  and  is  made  to  be  coincident 
with  the  diffracted  information  from  the  AOTDL  for  each  tap  position  and  for  a  given  input  frequency  (the 
AOTDL  center  frequency).  The  reference  beam  then  interferes  temporally  (heterodynes)  with  the  beam 
from  the  signal  path  to  form  the  desired  output  at  the  system  IF.  The  Doppler  shift,  due  to  the  AOSLM,  will 
cancel  out  when  the  two  paths  are  heterodyned. 

The  following  interpretations  of  the  math  will  be  useful  in  working  through  the  equations  to  follow: 

1 .  Tilted  plane  wave:  exp  [jox]  or  exp  [jat],  where  x  and  x  are  spatial  coordinates, 

2.  Magnified  and  shifted  image:  M‘^g(x/M  -  dx/M),  where  dx  is  the  positional  shift,  and 

3.  Time  delay  as  a  linear-in-frequency  phase  shift  for  a  tone: 
exp(j2rc/(t-T)J  =  exp[j2rt/tlexp[-j27i/T) . 

4.3.1 .1  AOSLM  Mathematical  Description 

We  first  consider  the  tap  function  input  to  the  AOSLM  given  by: 

N 

f(t)  =  2^ajCos2Tc^t  (4-19) 

i  =  1 

where  aj  and  f  \  are  the  tap  amplitude  and  frequency,  respectively.  The  term  aj^  has  units  of  power.  Note 
that  the  phase  of  the  AOSLM  input  will  not  affect  the  processor  output  since  it  will  cancel  during  heterodyne 
detection.  This  signal  results  in  the  acoustic  wave  diffracted  output  for  +1  diffraction  (ignoring  the 
undiffracted  beam  blocked  at  plane  PI)  given  by: 

N 

where  A  is  the  incident  plane-wave  amplitude,  r\i  is  the  AO  cell  diffraction  efficiency  in  percentage  per  RF 
watt,  and  w(x)  is  the  apodization  function  due  to  beam  shape,  acoustic  aperture,  acoustic  attenuation,  lens 
aberrations,  etc.  The  nonuniform  frequency  response  of  the  AOSLM  is  included  for  convenience  in  the 
coefficients  aj.  The  frequency  response  can  be  electronically  equalized  by  use  of  a  pre-multiplication  look¬ 
up  table  that  is  a  function  of  /, 
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As  described  above,  lenses  Li.  L2  and  L3  form  the  Fourier  transform  of  Equation  4-20  and 
display  this  on  the  aperture  of  the  AOTDL  witfi  iiequency  coordinates  given  by  Equation  4-17. 
Mathematically,  this  signal  is  given  at  Plane  Pi  as  (Reference  19): 


F(t.t)  = 


N 


(X)  ajexpfjajt^ft- 1  -  ^llexp 

i  =  i  L  ^ 


-'\2nx 


1J 


dx 


(4-21) 


where  the  k  phase  shift  (l/j)  is  neglected.  Using  substitutions  and  Fourier  transform  properties,  we  arrive 
at  the  transform: 


F(t,T) 


i]exp[i2^/.(t  -  1)]. 


(4-22) 


After  lenses  L2  and  L3,  Equation  4-22  can  be  modified  to  include  the  telescope  magnification  to  yield  the 
resultant  illumination  of  the  AOTDL  given  by: 


N 


F  (t.T)  =  A 


"I'a"' 


1/2 


Xa^W 
=  1 


Tf, 


/i 


^*1*3 


exp[i27tyj(t  -  ^)]- 


4.3.1 .2  Reference  Beam  Path 
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The  pattern  F(t,r)  passes  through  two  legs  of  a  Mach-Zehnder  interferometer,  one  of  which 
contains  the  AOTDL  and  one  of  which  will  be  referred  to  as  the  reference  path.  The  input  to  the  AOTDL 

consists  of  individual  beams  for  i=l ,  2 . N.  all  of  which  are  parallel  to  each  other  and  individually  diveiging. 

This  is  a  consequence  of  the  Fourier  transform  properties,  where  points  on  the  input  map  to  plane  waves 
on  the  output,  and  plane  waves  on  the  input  map  to  diverging  parallel  beams  on  the  output.  The  final  lens, 
L4,  then  transforms  the  AOTDL  diffracted  and  undiffracted  beams  onto  the  detector  to  form  a  magnified 
version  of  the  beams  at  plane  P2,  i.e.,  lenses  L3  and  L4  act  as  another  telescope  with  magnification.  This 
is  also  effectively  the  image  of  the  AOSLM  aperture,  consisting  of  overlapped  plane  wave  beams  at  various 
incident  angles  corresponding  to  the  various  AOSLM  input  frequencies.  Therefore,  the  reference  beam  at 
the  detector  is: 


fj. 


*2*4  J 


S  a,exp 


i  =  1 


-j2n/: 


fifoxn 


exp 


(4-24) 
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4.3.1 .3  Signal  Beam  Path:  Signal  Carrying  Beam 


Now  cxjnsider  the  signal  path  of  the  Mach-Zehnder  interferometer.  Consider  an  input  signal  to 
the  AOTDL  given  by: 

S  (t)  =  a  (t)  cos  [2n7  t  +  <t)  (t)  ]  (4-25) 

which  is  a  complex  modulation  of  a  carrier,  For  a  better  understanding  of  the  equations  to  follow,  this 
signal  is  represented: 

M 

s(t)  =  J]bjCos[2rt?:t  +  <{i.].  (4-26) 

j  =  i 

Referring  back  to  Equation  4-17,  the  AOTDL  is  tapped  at  positions: 

Xi  =  (4-27) 

resulting  in  the  diffracted  signal  outputs  (recall  that  propagation  is  in  the  -  x  direction): 


(4-28) 


The  next-to-last  term  contains  a  component  exp  [  \2k'F^N^,  which  is  the  resulting  angular  offset  of  the 
diffracted  beam  with  respect  to  the  undiffracted  beam.  The  offset  at  the  input  center  frequency  is  removed 
relative  to  the  reference  beam  of  Equation  4-24  by  tilting  the  mirror  in  the  signal  path  and  aligning  the  beam 
combiner  cube  so  that  the  signal  and  reference  paths  superimpose  after  recombination,  as  shown  in 
Figure  4-23.  We  then  define  the  signal  frequencies  relative  to  the  center  frequency  as  Ajj  (A9j=0  for  the 
center  frequency).  Therefore,  after  removing  this  offset. 


S(t,x) 


M 

i=  1 


^2 

—  biexp[j(^J 


(4-29) 


exp  i27t(?j+/j)  (t-I)]exp 

L  3 


Lens  L4  then  takes  the  Fourier  transform  of  this  signal  as  was  done  in  Equations  4-21  through 
4-23  (the  exponential  in  the  Fourier  integral  is  now  positive).  Therefore,  the  signal  beam  at  the  detector  is: 
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N  M 


I  I  "  ^jbjexp[j4,j] 

-  oo  i  —  1  j  —  1 


f\ 


'1‘3  *3 


•  exp  j2n(7j+/j)  (t-I)]exp  exp  dx. 

_  3_  _  4_ 


This  can  be  simplified  by  using  the  Fourier  transform  properties  to  yield; 


s(t,x)  =  A 


f.f,  ^  /xfif-i  fifo\  r 

rr'if  Z  Z^i'^lrr  pp 

V4  jt^  '^’2'4  ”2''a^  L  V4''aJ 


XAfrf,f3 


•  bjexp[j<|>.]exp  -\2nf^ - exp[j2n(?-j  +  y:)  (t- -)]. 

‘2''a 


To  Simplify  this  notation,  we  define  a  position  independent  phase  offset  given  as: 


A<I)..  =  -211/ 


XAjfifg 


(4-30) 


(4-31) 


(4-32) 


Which  is  dependent  on  both  tap  position  and  signal  frequency.  This  phase  term,  which  is  linear  in  tap 
position  and  linear  in  offset  signal  frequency  (relative  to  the  signal  center  frequency),  is  apparently  the 
relative  time  delay  of  the  signal  for  a  given  tap  position,  as  was  shown  in  Item  3  immediately  preceding 
Subsection  4.3.1 .1 .  This  apparent  time  delay  is  the  source  of  the  limitation  of  this  architecture,  as  will  be 
described  in  Subsection  4.3.1 .7.  The  signal  beam  at  the  detector  is: 


N  M 


i  =  i  j=  1 


’2’4''a 


(4-33) 


’  bjexp  [j  (<|>j  -t-  AO,.)  1  exp[^j2ji  (?-.  +/.)  (t  -  ^)]  • 


It  is  instmctive  at  this  point  to  compare  Equations  4-33  and  4-24.  We  see  that  these  two  detector 
illuminations: 


*4 

1 .  Differ  in  position  by  an  amount  XA  ?j  — , 

'^a 

2.  Differ  in  temporal  modulation  by  the  signal  modulation  given  in  Equation  4-26,  and 

3.  Have  a  tap  position  and  signal  frequency  offset  dependent  phase  shift  of  the  temporal 
rrwdulation. 
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Number  3  above  is  key  to  the  delay  line  properties  of  the  system.  Note  that  if  A?j  is  fixed,  the  phase  delay 
is  a  linear  function  offset  of  tap  position,  as  it  should  be.  Also,  a  constant  delay  will  result  in  a  phase  shift 
that  is  linear  with  signal  frequency.  Therefore,  it  appears  that  the  linear  phase  given  in  Equation  4-32 
corresponds  to  the  time  delay  of  the  signal.  This  does  not  actually  occur,  as  will  be  described  in  Subsection 
4.3.1 .7. 

4.3.1 .4  Signal  Beam  Path:  Undiffracted  Beam 


The  undiffracted  beam  from  the  AOTDL  will  be  identical  to  the  reference  beam  in  the  other  path 
except  for  the  positional  offset  introduced  by  the  tilted  mirror.  This  is  equivalent  to  the  signal  beam 
mathematics  with  A?j  replaced  by  9c-  The  undiffracted  beam  is  therefore  written  as: 


M 

'-nf  I'^f 

V  i=  1  ^ 


.1/2 

^  xf,f 
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p 

expljAO'.-iexp 
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where  AO'jj  is  given  by  Equation  4-33  with  ^9'^  replaced  by  9c 

4.3.1 .5  The  Detected  Heterodyned  Signal 


The  detected  signal  for  the  three  incident  amplitude  waves  is  given  by  the  square-law  summation: 
L/2 

d(t)  =  J  If  (t,x) +s(t,x) -)-u(t,x)|^dx  (4-35) 

-L/2 


where  L  is  the  size  of  the  detector,  and  the  three  signals  are  given  in  Equations  4-24,  4-33,  and  4-34.  This 
signal  can  be  expanded  as: 


L/2 

d(t)  =  j  [lf(t,x)|^ -t-|s(t,x)|^■^|u(t,x)|^-h2Re  {f*  (t,x)s(t,x) } 
-L/2 

-t-2Re  {f*  (t,x)u(t,x) }  -t-2Re  {s*  (t,x)u(t,x) }  ]dx. 


(4-36) 


The  fourth  term  contains  the  desired  output  information  at  the  IF,  9'c-  Note  that  the  first,  second,  third,  and 
fifth  terms  all  represent  biases  that  are  temporally  at  dc,  and  thus  can  be  filtered  out  after  detection  using  a 
BPF  centered  at  the  signal  IF.  The  form  of  the  desired  signal  output  is  first  derived  and  then  used  to  calculate 
the  form  of  the  sixth  term. 
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The  term  containing  the  desired  output  signal  is  written  as: 


L/2 


d^(t)  =  J  2Re{f*  (t.x)s(t.x)}dx 


-L/2 
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(4-37) 


•  exp  I  j  (Oj  +  A<t)^j) )  exp  [j27t  -f/^)  (t 


This  equation  consists  of  the  product  of  two  summations,  one  over  tap  position  i  and  the  other  over  both  tap 
position  k  and  signal  frequency  component  j.  This  will  result  in  products  for  which  i  in  the  first  sum  equals  i 
in  the  second  sum,  and  for  which  i  is  not  the  same  as  k  in  either  sum.  The  case  for  which  i=k  is  the  desired 
case,  whereas  the  cross-terms  result  in  undesired  harmonics  of  the  input  signal  offset  in  frequency. 

Consider  i  =  k  in  both  summations  written  as: 


‘2»4  '  J  ^  ^2UJ  ^  V4  ^‘2''a 
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1j  =  1 

■  ai^bjexp [j (<t).  +  AOj.) ) exp|^j27t!r. (t -  I)]dx 


(4-38) 


From  linearity,  the  integral  of  the  real  part  of  a  signal  is  equal  to  the  real  part  of  the  integral,  and  the  integral 
of  the  sum  is  equal  to  the  sum  of  the  integral.  Therefore  Equation  4-38  can  be  rewritten  to  yield  the  desired 
signal  output: 


d4(t;i=k)  ^  2a2|1^ti3/2 
'2’4 


(4-39) 
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For  each  tap  position  i,  the  appropriate  delayed  version  of  the  signal  (the  delay  is  due  to  the  phase  term 
A(l>ij),  can  be  reconstructed  multiplied  by  the  square  of  the  tap  amplitude,  and  weighted  with  the  integral 
shown,  which  is  dependent  on  the  signal  frequency. 

The  term  for  which  the  i  s  are  not  equal  is  given  by  the  general  form  in  Equation  4-37.  These 
cross-  terms  are  due  to  the  interference  of  one  tapped  signal  output  with  the  undiffracted  beam  for  a 
neighboring  tap.  For  this  case,  the  equation  can  be  written  to  express  the  undesired  cross  terms  as: 
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The  undesired  signal  output  is  at  a  carrier  frequency  offset  from  the  input  signal  carrier  frequency  by  an 
amount  equal  to  the  difference  of  the  tap  Doppler  offsets.  Therefore,  this  carrier  frequency  will  fall  within 
the  system  pass  band  if  the  taps  are  closely  spaced.  This  term  is  weighted  by  the  integral,  which  can  be 
viewed  as  the  integral  over  a  limited  detector  extent  of  the  product  of  offset  AOSLM  imaged  apertures. 
Assume  for  simplicity  that  the  product  of  the  weighting  functions  windowed  by  the  detector  aperture  is  a  rect 
function,  which  assumes  total  overlap  on  the  detector  which  is  most  likely  not  the  case.  Then  Equation  4- 
40  can  be  expressed  as  (Reference  20): 
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which  reduces  to 
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where  sirK;(P)  =  sinnp/np.  As  an  example  of  the  effect  of  this  equation,  the  frequency  separation  required 
to  reach  the  first  null  of  the  sine  function,  which  occurs  for  p=l,  is  calculated.  For  representative  focal 
lengths  (fi,f2.f3.f4)  =  (.5, .05,. 3, .3),  detector  size  L=1  mm.  and  Vg  =  4200m/s,  then/j  -/k=  420  kHz.  Observe 
that  the  larger  the  detector,  the  closer  the  tap  frequencies. 

These  cross-terms  are  due  to  the  interference  of  different  tap  positions  with  each  other,  and  are 
due  in  part  to  the  spatial  coherence  of  the  tap  input  (all  the  light  comes  from  a  single  laser).  CPI  Keefer’s 
investigation  of  laser  diode  array  inputs  to  the  AOTOL  would  be  a  potential  solution  to  this  problem.  If  the 
independently  modulated  lasers  in  this  array  are  spatially  incoherent,  then  adjacent  taps  will  not  interfere 
with  each  other.  The  possibility  of  driving  such  a  laser  diode  array  with  the  output  of  a  photodiode  array/ 
integrator  would  potentially  be  a  very  attractive  alternative  to  the  AOSLM  approach  when  considering  the 
cross-terms. 

Now  consider  term  six  in  Equation  3-36  due  to  the  interference  of  the  undiffracted  and  diffracted 
beams  from  the  signal  path.  This  term  can  be  written  as: 
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Once  again,  there  are  two  cases  depending  on  whether  i  =  k  or  i  ^  k.  Consider  only  the  general  case,  which 
includes  the  special  case  of  i=k,  for  which: 
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For  the  system  parameters  of  our  hardware  implementation,  the  physical  separation  between  the  two 
weighting  functions  in  the  integral  is  larger  than  their  width,  i.e.,  the  two  terms  appear  separated  at  the 
detector.  As  an  example,  for  the  focal  lengths  and  velocity  above,  (A;7j.?c)  =  (5MHz,40MHz)  and  a 
wavelength  of  532  nm,  the  separation  is  given  as  13.3  mm.  Thus,  no  interference  occurs  for  this  term. 

In  summary,  the  only  terms  that  impact  the  output  of  the  tapped  delay  line  filter  are  the  signal  term 
and  cross-terms  in  Equations  4-39  and  4-40.  The  detected  cross-terms  will  be  suppressed  as  the  tap 
frequency  separation  increases  and  as  the  detector  size  is  reduced. 

4.3.1 .6  Interferometric  Architecture  Experimental  Results 

A  Single-channel  AOTDL  architecture  was  fabricated  according  to  the  theory  described  above. 
This  system  operated  with  the  existing  AO  cells  and  at  a  40-MHz  IF.  Initial  detector  outputs  were  very  noisy, 
and  had  low-frequency  modulation  imposed  on  the  desired  output  at  40-MHz.  This  low-frequency 
modulation  was  due  to  the  interference  of  the  high-frequency  terms  with  each  other  to  induce  baseband 
modulations.  A  40-MHz  BPF  with  a  bandwidth  of  3  MHz  was  inserted  after  amplification  to  greatly  improve 
the  output  signal.  There  was  no  noticeable  difference  when  this  filter  was  placed  before  the  amplifier. 

Previous  implementations  were  not  set  up  to  achieve  the  imaging  and  Fourier  transform 
operations  described  above  for  lenses  Li  through  L4,  although  they  approximated  this  condition.  The 
system  was  carefully  realigned  to  achieve  this  operation,  with  new  lenses  being  inserted  where  appropriate. 
A  1-in.  aperture  collimation  tester  was  employed  to  accurately  collimate  the  system.  It  was  recommended 
that  a  cemented  doublet  be  employed  after  the  pinhole  in  the  initial  filtering  and  beam  expanding  stage  to 
reduce  aberrations.  This  greatly  improved  the  collimation  performance.  All  plano-convex  lenses  were 
oriented  so  their  planar  side  faced  the  converging  or  diverging  part  of  the  beam.  Under  these  conditions, 
an  image  of  the  AOSLM  was  formed  on  the  detector.  The  size  of  this  image  was  changed  by  adjusting  the 
aperture  at  the  AOSLM  input.  A  small  spot  was  formed  in  the  AOTDL,  although  this  spot  was  often  distorted, 
most  likely  due  to  aberrations  primarily  induced  by  the  short  focal  length  lens  L2.  No  attempt  was  made  to 
modify  the  apodization  function  of  the  AOSLM  as  described  above  in  the  mathematical  development.  This 
apodization  would  result  in  smaller  tap  spots  in  the  AOTDL  and  generally  improve  the  performance  of  the 
system. 

The  alignment  of  the  Mach-Zehnder  interferometer  was  critical.  Initially,  the  noncommon  path 
included  the  lens  and  spatial  filter  slit  after  the  AOTDL,  used  previously  to  remove  the  undiffracted  beam, 
and  lens  L3.  These  components  were  set  outside  the  interferometer,  and  it  was  realized  at  this  time  that 
filtering  the  undiffracted  beam  was  unnecessary.  Every  attempt  was  made  to  reduce  the  size  of  the 
interferometer,  greatly  enhancing  the  stability,  which  was  quite  good  for  a  system  employing  tilt  tables  and 
adjustable  mirrors.  However,  during  testing  of  the  system’s  operation,  it  was  discovered  that  the  phase  of 
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the  output  signal  did  not  follow  the  tap  position  (as  selected  by  the  AOTDL)  as  expected.  Although  many 
modifications  were  attempted  in  the  geometry  of  the  interferometer,  the  desired  phase  behavior  was  not 
achieved.  We  suggested  the  addition  of  a  holographic  grating  in  the  same  location  as  the  AOTDL  but  in  the 
reference  path,  allowing  a  symmetric  geometry  to  be  obtained  and  eliminating  the  change  in  path  length 
difference.  Two  holographic  gratings  were  made  at  the  Photonics  Center  by  Stephen  Kupiak  of  the 
University  of  Alabama  in  Huntsville,  but  efficiency  was  not  very  high  and  time  did  not  permit  a  full  evaluation 
of  the  system  after  insertion  of  these  gratings.  Subsection  4.3.1 .7  develops  the  source  of  this  anomalous 
phase  behavior,  and  a  modified  architecture  that  overcomes  this  limitation  is  the  subject  of  Subsection  4.3.2. 

The  output  at  the  40-MHz  IF  was  qualitatively  compared  to  the  AOTDL  input,  and  it  was  found 
that  there  was  very  little  distortion  of  the  input  signal.  Figure  4-24  shows  such  an  output  for  a  DSB-SC 
modulation  of  the  40-MHz  carrier  with  a  l-MHz  tone.  Although  this  plot,  created  on  1  August  1991 ,  shows 
a  2-mV/div  scale  for  the  optical  signal,  outputs  of  greater  than  100  mV  were  typically  achieved  in  the 
following  weeks.  The  canceled  signal  is  also  shown  in  this  figure.  The  cancellation  of  the  signal  at  the  IF 
is  critical  to  meeting  the  needs  of  the  radar  system,  but  this  requirement  introduces  other  real-world 
complications  not  present  with  the  previous  base-band  implementations.  Because  the  signals  are  being 
subtracted  on  the  carrier,  the  carriers  must  be  aligned  in  phase,  in  addition  to  the  alignment  of  the  envelope 
modulations.  If  the  tap  position  is  moved  to  achieve  a  180°  phase  shift  of  the  carrier,  no  noticeable  change 
in  the  position  of  the  envelope  will  result,  but  the  cancellation  will  go  from  good  cancellation  depth  to 
constructive  interference  (a  doubling  of  the  signal  amplitude).  Therefore,  it  is  clear  that  precise  control  of 
the  tap  position  is  required.  It  also  becomes  clear  why  the  phase  is  important  in  the  cancellation.  If  the 
auxiliary  channel  of  a  single-channel  system  is  phase  shifted  (a  phase  shift  of  the  carrier  but  no  time  shift 
of  the  envelope)  relative  to  the  main  channel,  then  to  get  the  best  cancellation,  the  envelopes  must  be 
aligned  and  the  phase  shift  must  be  compensated. 

Although  the  frequency  response  (magnitude  and  phase)  of  the  system  was  not  fully 
characterized,  the  magnitude  of  the  output  v  as  measured  as  a  function  of  tap  position  (as  selected  by  the 
AOSLM)  for  a  given  alignment  condition  and  a  single-tone  input  to  the  AOTDL.  The  results  of  this  are  shown 
in  Table  4-3,  where  it  is  seen  that  the  magnitude  response  is  very  flat  over  the  span  of  AOSLM  input 
frequencies.  This  is  especially  surprising  considering  the  composite  effects  of  the  AOSLM  frequency 
response,  the  AOTDL  attenuation  properties,  and  the  heterodyne  efficiency  of  the  process.  Although  this 
result  is  very  promising,  this  measurement  was  made  for  one  condition  of  system  alignment,  and  may  be 
different  when  the  system  is  optimally  aligned  for  magnitude  and  phase  response. 

Finally,  the  generation  of  the  undesired  cross-terms  (described  by  Equation  4-40)  was 
qualitatively  examined  and  found  to  adhere  approximately  to  theory.  The  AOSLM  input  consisted  of  a  DSB- 
SC  modulation  of  a  40-MHz  carrier  with  a  tone.  This  tone  was  decreased  in  frequency  from  greater  than 
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Figure  4-24.  Representative  Tapped  Delay  Line  Signals 
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Table  4-3.  Amplitude  Response  Versus  Tap  Position 


5  MHz  as  the  detector  output  was  examined  on  the  spectrum  analyzer.  As  the  tone  approached  0.5  to  1 .0 
MHz,  terms  at  the  difference  frequencies  7^  +/k-/i  Equation  4-40  were  apparent.  These  terms  grew  in 
magnitude  as  the  tone  was  further  decreased,  as  expected  No  characterization  of  the  effect  of  detector 
size  on  cross-term  suppression  was  performed. 

4.3.1 .7  Interferometric  Architecture  Limitations 

When  this  Mach-Zehnder  architecture  was  developed,  it  was  experimentally  observed  that  it  was 
very  difficult  to  make  the  signal  information  delay  the  appropriate  amount  when  the  tap  position  was 
adjusted  by  tuning  the  input  to  the  AOSLM.  In  fact,  it  appeared  that  when  the  system  was  in  its  best  state 
of  alignment,  the  carrier  was  stationary  even  when  the  modulation  envelope  moved  at  the  appropriate  rate 
relative  to  the  reference  signal  (RF  input  to  the  AOTDL).  When  the  AOTDL  was  moved  relative  to  the  tap 
position,  true  time  delay  was  consistently  observed.  This  implied  that  moving  the  AOTDL  relative  to  the  tap 
was  significantly  different  from  moving  the  tap  relative  to  the  AOTDL. 

The  arrangement  shown  in  Figure  4-25  is  first  considered.  In  this  figure,  a  single  frequency  is 
input  to  the  AOTDL  and  a  plane  wave  illumination  is  provided.  A  reference  beam  is  made  to  be  coincident 
with  the  diffracted  beam,  in  a  similar  fashion  to  the  reference  beam  shown  in  Figure  4-23  The  diffracted 
beam  and  a  small  reference  beam  in  this  case  are  collimated  and  colinear,  therefore  having  a  constant 
phase  front  across  the  overlapping  beams.  Placing  a  small  detector  in  the  beam  and  moving  this  detector 
orthogonal  to  the  propagation  direction  (while  at  the  same  time  moving  the  reference  beam  to  overlap  on 
the  detector)  results  in  no  relative  position-dependent  phase  shift  of  the  heterodyned  signal.  This  is 
equivalent  to  tapping  the  AOTDL  at  different  positions  and  observing  the  output  of  the  heterodyne  process 
relative  to  the  RF  input  to  the  AOTDL. 

Next,  an  identical  geometry  was  considered  and  the  AOTDL  was  moved,  but  the  detector  position 
and  reference  beam  position  were  fixed.  As  for  a  moving  grating,  the  phase  of  the  output  will  change  relative 
to  the  input  as  the  AOTDL  is  moved.  One  conceptual  way  of  observing  this  phase  shift  is  to  consider  the 
imaging  of  a  grating.  As  the  grating  is  moved,  the  image  will  rrwve;  therefore,  a  small  detector  will  see  the 
light  go  from  dark  to  bright  as  the  fringes  pass  by.  For  a  heterodyne  system,  this  is  equivalent  to  seeing  the 
phase  of  the  heterodyned  output  change  linearly  with  AOTDL  position. 

The  case  of  concern  (moving  tap  position)  can  be  mathematically  described  with  a  simple 
example  of  a  pulsed  carrier  The  pulse  input  is  described  as: 

s(t)  =  p^(t)exp[j2nyi  (4-45) 

where  x  denotes  the  pulsewidth  and  f<;  is  the  carrier  frequency.  Assuming  that  t=0  corresponds  to  the  centei 
of  the  AOTDL  (x=0),  and  T  is  the  AOTDL  time  aperture,  the  pulse  can  be  written  as: 
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(4-46) 


X 


Adding  the  relerence  beam  that  travels  in  the  same  direction  as  the  center  frequency  to  this  signal,  and 
heterodyne  detecting,  an  output  signal  was  obtained,  given  as  a  function  of  time  and  tap  position  by: 
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Now.  for  x=0  (a  tap  position  at  the  center  of  the  AOTDL),  the  output  is  given  by: 


l(t,0)  =  biases  +  2Re|  P^(t-^)exp[^j2rtfj.(t-^)]i 

(4-48) 

and  for  x=  (note  that  L=VaT),  the  output  is  given  by: 

l(t,I)  =  biases  +  2Re|  p^(t)exp[j2icfg(t  -  ^)]  . 

^  (4-49) 

Therefore,  two  tap  positions  see  the  envelope  of  the  waveform  at  the  different  times  but  the  carrier  phase 
relationship  to  the  Dulse  is  not  constant.  This  resolves  the  apparent  paradox  observed  in  the  hardware 
implementation,  namely  that  the  envelope  of  the  signal  moved  appropriately  but  the  carrier  moved  relative 
to  the  envelope  and  even  appeared  stationary. 

It  is  interesting  at  this  point  to  refer  back  to  the  mathematics  of  the  Mach-Zehnder  architecture. 
A  key  term  to  be  analyzed  is  the  tap  position  and  AOTDL  frequency-dependent  phase  shift  of  the 
heterodyned  signal.  This  is  defined  in  Equation  4-32  as: 

A.AT'.fJ,, 

AO..  =  -271/ - (4-50) 

'I  f  v2 

where /i  and  A?j  are  the  AOSLM  input  frequency  (defining  the  tap  position)  and  frequency  offsef  of  the  signal 
from  the  AOTDL  center  frequency.  Note  that  the  AOTDL  center  frequency  is  defined  for  this  case  as  the 
frequency  that  exactly  aligns  with  the  reference  beam  of  the  Mach-Zehnder  interferometer.  This  was  initially 
interpreted  as  the  desired  phase  shift,  which  is  linear  in  frequency  (thought  to  result  in  true  time  delay)  and 
linear  in  tap  position.  The  key  problem  with  this  interpretation  is  that  the  phase  shift  is  a  function  of  the 
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difference  frequency.  A?],  and  not  the  AOTDL  input  frequency.  For  an  input  equal  to  the  center  frequency, 
the  phase  shift  is  completely  independent  of  tap  position  (note  the  similarity  to  the  argument  in  reference  to 
Figure  4-25).  Equation  4-50  above  preserves  the  true  time-delay  of  the  envelope  of  the  modulation  but 
does  not  time-delay  the  carrier. 

A  straightforward  way  to  mathematically  understand  this  phenomenon  is  to  consider  a  signal  on 

a  carrier. 


s(t)  =  a(t)exp[j27ty]  (4-51) 

which  has  a  Fourier  transform: 

S(f)  =  A(f-f^).  (4-52) 

A  true  time-delay  of  s(t)  is  written  as. 

s(t-T)  =  a(t-x)exp[j2nf^(t-T))  (4-53) 

which  has  a  Fourier  transform; 

S(f,T)  =  exp[-j2nxf]  A(f-f^)  (4-54) 

By  imposing  a  linear  phase  as  a  function  of  difference  of  frequency,  as  in  Equation  4-50,  it  can  now  be 
written; 

S'(f.x)  =  exp[-j2jix(f-f^)l  A(f-f^)  (4-55) 

which  is  in  the  time  domain  given  by: 

s'(t-x)  =  exp[j27txf^]a(t-x)exp[j27tf^(t-x)]  =  a(t-x)exp[j2ny] .  (4-56) 

The  carrier  for  this  case  is  seen  to  be  stationary  with  a  time-shifted  envelope. 


To  remedy  the  situation,  it  can  be  observed  that  to  convert  to  the  AOTDL  input  frequency  a 
reference  beam  must  be  employed  that  appears  to  emanate  from  the  tap  at  an  angle  equivalent  to  the 
undiffracted  beam.  This  is  equivalent  to  shifting  the  position  of  the  reference  beam  to  align  with  the 
undiffracted  beam  rather  than  to  align  with  the  center  frequency  beam.  Because  of  this  realization,  it  is 
apparent  that  the  reference  beam  is  no  longer  needed  and  the  undiffracted  beam  can  be  used  instead.  This 
forms  the  basis  of  the  noninterferometric  architecture  described  in  Subsection  4.3.2. 
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4.3.2  Noninterferometric  Architecture  Design  and  implementation 


Several  practical  issues  relate  to  the  use  of  the  undiffracted  beam  as  a  reference  beam  in  the 
AOTDL  filter.  First,  the  degree  of  overlap  of  the  two  beams  at  the  detector  is  very  small  since  the  angles  at 
the  AOTDL  are  very  large.  To  remedy  this  situation,  the  tap  must  be  made  very  small,  resulting  in  a  larger 
beam  on  the  detector  for  better  overlap  (refer  back  to  Equation  4-44).  Secondly,  the  amplitudes  of  the 
diffracted  and  undiffracted  beams  are  coupled,  thus  the  process  is  not  linear.  A  technique  for  achieving  both 
better  overlap  and  less  coupling  of  power  between  the  reference  and  signal  beams  is  to  use  an  electronic 
reference.  This  architecture,  which  was  designed,  fabricated,  and  tested  during  this  effort  will  now  be 
described.  The  resulting  architecture  is  similar  in  layout  to  prior  architectures  reported  in  Reference  1 . 

The  noninterferometric  tapped  delay  line  architecture  is  shown  in  Figure  4-26.  The  system  is 
equivalent  to  the  interferometric  approach  shown  in  Figure  4-23  with  the  reference  path  removed  and  an 
electronic  reference  added.  Observe  that  this  architecture  is  effectively  an  AO  point  modulator  with  a 
moving  tap  position. 

The  mathematics  for  the  noninterferometric  architecture  closely  parallels  that  of  the 
interferometric  architecture. 

4.3.2.1  SignaUCarrying  Beam  and  Undiffracted  Beam 


For  the  noninterferometric  architecture,  the  signal-carrying  beam  and  the  undiffracted  beam  are 
identical  to  those  for  the  interferometric  case  except  that  the  frequencies  are  relative  to  dc  and  not  7^. 
Therefore,  the  signal-carrying  beam  can  be  written  as  (see  Equation  4-33): 
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where  AOjj  is  now  given  as; 


AO..  =  -2nf 

I)  -'i 
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(4-58) 


Likewise,  the  undiffracted  beam  can  be  written  as; 
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(4-59) 
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where  AO'jj  =^^0=  0  in  Equation  4-34  and  a  will  be  defined  later.  Note  that  this  beam  is  alnxjst  identical  to 
the  reference  path  beam  given  by  Equation  4-24,  as  it  should  be. 

4.3.2.2  Addition  of  a  Diffracted  Electronic  Reference  Beam 


For  reasons  to  be  explained  below,  related  to  the  position  offset  of  the  diffracted  and  undiffracted 
beams  at  the  photodetector,  an  electronic  reference  is  ar>plied  and  given  by: 

r(t)  =  acosl2jtfrg^t]  (4-60) 

This  electronic  reference  beam  behaves  as  the  diffracted  signal  beam  given  in  Equation  4-34;  therefore,  the 
electronic  reference  beam  can  be  written  as; 


+  exp  -j27tf 
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•  aexp[jA<t).g^]exp[j2jr(?-g^  +  /j)  0-^)] 

where  A<I»ier  is  given  by  Equation  4-58  with  replaced  by  fer 
4.3.2.3  The  Detected  Heterodyned  Signal 


(4-61) 


The  detected  signal  for  the  three  incident  amplitude  waves  is  given  by  the  square-law  summation; 
L/2 

d(t)  =  I  |u(t,x) +s(t,x) -i-r(t,x)l^dx  (4-62) 

•\J2 

where  L  is  the  size  of  the  detector  and  the  three  signals  are  given  in  Equations  4-57, 4-59,  and  4-61 .  This 
signal  can  be  expanded  as  . 


4-66 


(4-63) 


L/2 

d(t)  =  J  [|u(t,x)|^  +  |s(t,x)|^  +  lr(t,x)|^  +  2Re{u*  (t.x)s(t,x)} 

-\J2 

+  2Re  {u*  (t.x)  r(t.x) }  +  2Re  {s*  (t.x)  r(t,x) }  Idx. 

The  fourth  term  comes  out  of  the  detector  after  heterodyning  at  the  center  frequency,  %  whereas  the  sixth 
term  is  at  Ter-  Both  terms  contain  the  desired  output,  but  to  generate  the  output  at  the  processor  IF, 
the  sixth  term  must  be  mixed  to  the  IF  with  the  frequency  Note  that  the  first,  second,  and  third  terms  all 
represent  biases  that  are  temporally  at  dc,  and  can  be  filtered  out  after  detection  using  a  BPF.  Likewise, 
the  carrier  that  is  generated  by  the  fifth  term  at  can  be  removed  through  filtering.  The  calculated  form  of 
the  desired  signal  output  for  both  the  fourth  and  sixth  terms  is  given  below.  The  term  containing  the  desired 
output  signal  at  the  processor  IF  is  written  as; 

L/2 

d^(t)  =  j  2Re{u*  (t,x)s(t,x)}dx 

■U2 


•  exp 


(4-64) 


This  equation  consists  of  the  product  of  two  summations,  one  over  tap  position  i  and  the  other  over  both  a 
tap  position  k  and  signal  frequency  component  j.  This  results  in  products  for  which  i  in  the  first  sum  equals 
k  in  the  second  sum,  and  for  which  i  is  not  the  same  in  either  sum.  The  case  of  i=k  is  the  desired  case, 
whereas  the  cross-terms  result  in  undesired  harmonics  of  the  input  signal  are  offset  in  frequency. 


First  consider  i=k  in  both  summations: 
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From  linearity,  the  integral  of  the  real  part  of  a  signal  is  equal  to  the  real  part  of  the  integral,  and  the  integral 
of  the  sum  is  equal  to  the  sum  of  the  integral.  Therefore,  Equation  4-65  can  be  rewritten  to  yield  the  desired 
signal  output  at  the  processor  IF; 
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It  can  be  seen  that  for  each  tap  position  i,  the  appropriate  delayed  version  of  the  signal  has  been 
reconstructed,  (the  delay  is  due  to  the  phase  term  A4>,j),  multiplied  by  the  square  of  the  tap  amplitude,  and 
weighted  with  the  integral  shown,  which  is  dependent  on  the  signal  frequency. 

The  term  for  which  the  i  is  not  equal  to  k  is  given  by  the  general  form  in  Equation  4-64.  These 
cross-terms  are  due  to  the  interference  of  one  tapped  signal  output  with  the  undiffracted  beam  for  a 
neighboring  tap.  For  this  case,  this  equation  can  be  rewritten  to  express  the  undesired  cross-terms  as; 
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The  undesired  signal  output  is  at  a  carrier  frequency  offset  from  the  input  signal  carrier  frequency 
by  an  amount  equal  to  the  difference  of  the  tap  Doppler  offsets.  Therefore,  carrier  frequency  will  fall  within 
the  system  pass  band  if  the  taps  are  closely  spaced  in  frequency  input  to  the  AOSLM.  This  term  is  weighted 
by  the  integral,  which  can  be  viewed  as  the  integral  over  a  limited  detector  extent  of  the  product  of  AOSLM- 
imaged  AOSLM  apertures. 

The  sixth  term  in  Equation  4-63,  is  due  to  the  interference  of  the  electronic  reference  beam  and 
diffracted  signal  beam,  and  can  be  written  as: 


L/2  L/2 

dgft)  =  J  2Re{s*  (t.x)r(t,x)  }dx  =  j  2Re{s(t,x)r*  (t.x)  }dx 

-L/2  -L/2 


^  ^  /Xf,f3 


-L/2  U=1|=1 


*1*3 


>2'4 


'2''a 


-j2n/ 


*1*3X1 

*2*4''aJ 


bjexp[j(0j  +  A*|)lexp[j27i(7|+/;)  ^ 

i2"/Kg[^«exp[-jA<b,^j9xp[-j2Tt(ir^^  (*- J)]  jdx- 


(4-68) 


Once  again,  There  are  two  cases  depending  on  whether  i  =  k  or  not.  For  i=k: 
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Differences  between  this  term  and  Equation  4-65  are  the  more  substantial  overlap  of  the  two  window 
functions,  providing  more  heterodyne  efficiency  (motivation  for  using  an  electronic  reference),  a  carrier 
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frequency  now  equal  to  rather  than  the  system  IF.  and  an  additional  linear  phase  term,  A<I>jer  To 

obtain  this  desired  term  at  the  IF,  mixing  is  done  after  detection  with  the  electronic  reference  signal. 

Next  consider  the  general  case  of  i^^k  and  write; 
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These  terms  result  in  undesired  cross-terms  as  for  the  interferometric  case. 

4.3.2.4  Experimental  Results  and  Characterization 

On-site  at  the  Photonics  Laboratory,  the  interferometric  architecture  was  modified  to  create  the 
above-described  noninterferometric  architecture,  and  the  electronic  reference  was  inserted.  This 
subsection  documents  the  hardware  measurements  made.  Figures  4-27 and  4-28  st)ov/  the  optical  and 
electronic  systems  as  configured  in  the  laboratory  for  2-channel  operaiton. 

In  addition  to  configuring  the  system  for  noninterferometric  operation,  multichannel  AO  cells  were 
integrated  into  the  system  and  the  appropriate  lenses  were  inserted  betv  een  the  AOSLM  and  the  AOTDL. 
The  use  of  multichannel  devices  requires  that  the  output  of  the  AOSLM  be  transformed,  as  before,  in  the 
direction  of  acoustic  propagation,  but  in  addition,  the  orthogonal  direction  must  be  imaged  to  provide 
multiple  SLM  channel  illumination  of  the  multiple  channels  of  the  AOTDL.  A  number  of  configurations  were 
attempted  for  this  case  and  the  most  desired  approach  will  now  be  described.  With  equal  spacing  of 
acoustic  channels  for  the  two  AO  cells  selected,  the  imaging  must  be  1  ;1  in  the  direction  orthogonal  to  the 
acoustic  propagation  direction.  For  the  direction  of  acoustic  propagation,  the  system  must  Fourier  transform 
the  AOSLM  output  and  display  this  across  a  S-ps  aperture  at  the  AOTDL.  For  the  velocity  of  the  AOTDL, 
Va=4200  m/s,  this  corresponds  to  an  aperture.  Ax,  of  2.1  cm.  The  spread  in  frequencies  in  the  AOSLM  is 
nominally  set  at  40  MHz,  resulting  in  a  spatial  frequency  difference,  A/,,  of  9523.8  cycles/m.  Using  the 
equation. 

AX  =  Af^Xf  (4-71) 
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Figure  4-27.  Optical  Layout  for  Multichannel  AOTDL  Architecture 
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Figure  4-28.  RF  Layout  for  Multichannel  AOTDL  Architecture 
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it  is  determined  that  the  focal  length  requirement  is  f=4.145  m.  An  intermediate  Fourier  transform  can  be 
formed  arxJ  imaged  this  with  magnification  to  satisfy  the  equation  f  =  4.145  m.  As  an  example,  fi  and  fa 
can  be  made  to  equal  300  mm,  with  the  cylindrical  lens  having  a  focal  length,  fy  =  37.76  mm,  to  achieve  the 
desired  results.  For  the  system  fabricated  in  the  laboratory,  this  condition  was  approximated  by  having  f^ 
and  fa  equal  to  300  mm,  with  the  cylindrical  lens  having  a  focal  length,  fy  =  35.2  mm. 

Three  examples  of  spectrum  analyzer  outputs  for  the  system  are  shown  in  Figures  4-29  through 
4-31.  In  each  of  these  figures,  the  RF  input  into  the  AOTDL  and  the  optical  output  due  to  the  beating  of  the 
electronic  reference  beam  with  the  signal  beam  and  subsequent  mixing  back  to  the  system  IF  are  shown 
with  the  cancelled  output.  The  three  input  modulations  consisted  of  a  single  tone,  a  two-tone  waveform  with 
separation  of  3.2  MHz,  and  one  with  separation  10.0  MHz. 

Figure  4-32  (a)through  (f)  shows  the  spectrum  analyzer  display  of  the  AOTDL  input  (a  tone  at 
77  MHz  with  an  electronic  reference  at  63.6  MHz)  and  the  output  of  the  system  for  AOSLM  tap  frequencies 
of  60,  70,  79,  90,  and  100  MHz.  A  second  example  is  shown  in  Figure  4-33  (a)  through  (c).  It  was  of 
interest  to  trace  the  origin  of  the  different  spurious  tones  present  in  Figure  4-33(b}  to  isolate  them  and 
reduce  their  impact. 

The  output  of  the  desired  term  from  the  photodiode  is  given  by  the  difference  frequency,  fAoiDL 
-  fER.  for  a  signal  input  consisting  of  a  single  tone,  fAOTDL-  ^^6  output  of  the  mixer,  having  fER  as  its  input, 
is  given  by  two  tones  at  fAOTOL  desired  output)  and  2  fER  -  Iaotpl-  was  also  found  that  a  portion  of 
the  Iaotol  signal  reached  the  mixer  through  the  two  Pasternak  PE  2000  splitters  shown  in  Figure  4-33(c). 
This  resulted  in  the  generation  of  mixer  outputs  at  fER  and  2fAOTDL  '^er  aod  also  results  in  the  undesirable 
cross-talk  that  enters  the  photodetector  output  when  no  light  is  incident.  All  of  these  spurious  signals  have 
been  minimized  by  careful  shielding  and  the  addition  of  a  filter  that  blocks  the  feedthrough  of  the  fAOTOL 
signal  to  the  mixer.  The  final  measurements  that  were  made  utilized  an  HP  network  analyzer  configured  for 
8^2  paraineter  measurement  as  shown  in  Figure  4-34.  This  source  was  swept  across  a  wide  bandwidth  to 
measure  the  frequency  response  of  the  system  across  the  75  to  85  MHz  pass  band.  Two  examples  of  the 
response  are  shown  in  Figures  4-^  and  4-36  tor  electronic  reference  frequencies  of  69.4  and  7.28  MHz, 
respectively.  It  is  observed  in  both  cases  that  the  response  trails  off  as  the  AOTDL  input  separates  further 
from  the  electronic  reference.  This  is  due  to  the  overlap  of  the  two  beams  on  the  photodetector  as  given  in 
Equation  4-69.  This  frequency  response  can  be  improved  by  forming  better  spots  on  the  AOTDL  input 
(thereby  increasing  the  AOSLM  image  size  at  the  photodetector)  and  by  using  a  compensating  electronic 
filter  after  the  mixer.  The  improvement  of  the  optical  system  design  should  significantly  level  off  at  this 
frequency  response.  The  phase  response  for  a  representative  scenario  is  shown  in  Figure  4-37. 

Finally,  it  v/as  desired  to  measure  the  true  time-delay  properties  of  the  noninterferometric  tapped 
delay  line  architecture.  For  this  measurement,  the  output  of  the  system  after  mixing  up  to  system  IF  was 
examined  on  an  oscilloscope.  A  two-tone  waveform  with  a  4-MHz  tone  separation  on  an  80-MHz  carrier 
was  used  as  the  test  waveform.  The  test  waveform  was  used  as  a  trigger  reference  and  was  input  to  the 
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Figure  4-29.  Cancellation  of  a  Single  Tone 
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(b)  Optical  Output  for  60-MHz  AOSLM  Input 


Figure  4-32.  AODTL  System  Test  Measurements  (Case  1)  (Continued) 
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(d)  Optical  Output  for  79-MHz  AOSLM  Input 


Figure  4-32.  AODTL  System  Test  Measurements  (Case  1)  (Continued) 
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(f)  Optical  Output  for  100-MHz  AOSLM  Input 
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Figure  4-32.  AODTL  System  Test  Mecuurements  (Case  1)  (Concluded) 
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Figure  4-33.  AODTL  System  Test  Measurements  (Case  2) 
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(b)  Optical  Output  for  79-MHz  AOSLM  Input 
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Figure  4-33.  AOTDL  System  Test  Measurements  (Case  2)  (Continued) 
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Figure  4~34.  Network  Analyzer  Measurement  Setup 
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4-35.  Tapped  Delay  Line  Frequency  Response  for  f  =  69.4  MHz 


5dB/REF  -51.95  dB 
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igure  4-37.  Tapped  Delay  Line  Phase  Response  for  f  gp  ^  69.4  MHz 


oscilloscope's  channel  A  while  the  system  output  was  displayed  on  channel  B.  The  AOSLM  input  was  then 
tuned,  thereby  tuning  the  delay,  and  the  output  signal  was  observed  to  sweep  by  the  trigger  reference.  For 
this  waveform,  20  cycles  of  the  carrier  appear  between  nulls  of  the  envelope.  A  delay  sweep  from  null  to 
null  was  generated  and  the  carrier  was  seen  to  go  through  20  cycles,  thereby  verifying  the  true  time-delay 
behavior  of  the  system.  Further  examination  of  Equation  4-69  is  required  to  determine  the  effect  of  the  term 
-  AOjgp  and  how  true  time-delay  is  achieved.  This  further  analysis  should  also  include  the  post¬ 
detection  mixing  operation  to  achieve  the  output  at  the  system  IF. 

4.3.3  AOSLM  Intermodulation  Products 


A  significant  effect  in  this  tapped  delay  line  architecture  is  the  generation  of  spurious  signals  from 
the  AOSLM  due  to  diffraction  nonlinearities,  known  as  intermodulation  products.  For  two  input  frequencies, 
f  1  and  I2,  the  diffracted  output  will  correspond  to  these  two  frequencies,  but  two  outputs  in  the  pass  band  of 
the  device  at  frequencies  2fj-f2  and  2f2-fi  are  also  present.  This  is  shown  schematically  in  Figure  4-38. 
These  two  spurious  two-tone  intermodulation  products.  Ig,  are  suppressed  relative  to  the  intensity  of  each 
desired  frequency,  li  and  I2,  by  an  amount  equal  to; 
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For  !•)  =  Ig  =  Idesired'  f^ian. 


desired 
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desired 


where  rif  is  the  diffraction  efficiency  of  the  AOSLM  for  the  given  tap  frequency  power.  This  ratio  is  known 
as  the  spur-free  dynamic  range  of  the  system  and  defines  the  level  of  the  spurious  tones  relative  to  the 
maximum  two-tone  input  drive  powers.  This  can  be  thought  of  as  a  crosstalk  term  or  as  a  limitation  on 
contrast  ratio  when  multiple  taps  are  selected.  It  the  taps  are  normalized  so  that  the  sum  of  all  the  aj^  is 
equal  to  one,  then  this  product  is  included  in  the  equations  developed  above. 

This  can  be  used  in  the  design  to  set  a  maximum  drive  power  for  adequate  suppression  of  the 
two-tone  intermodulation  products.  This  drive  power  is  such  that  the  diffraction  efficiency  will  be; 


q,  =  6 


desired 


(4-74) 


For  example,  if  the  spurious-free  dynamic  range  of  the  system  is  required  to  be  40  dB,  then  Ig  /  Idesired 
equal  0.0001,  which  results  in  a  maximum  diffraction  efficiency  of  0.06,  or  6%  per  tap  frequency.  This 
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impacts  the  amplitude  of  the  output  signal  since  a  large  amount  of  the  light  into  the  system  is  lost  when  the 
undiffracted  AOSLM  output  is  filtered.  This  example  is  shown  on  the  plot  of  Figure  4-39,  which  is  for 
longitudinal  Te02.  A  0.06  light  output  level  is  squared  and  divided  by  36  to  find  the  two-tone  intermodulation 
product  which  will  be  at  6x10  ®,  as  shown  in  the  plot. 


SOURCE:  Reference  21 

TR-92-SI01-0027 

Figure  4-38.  Two-Tone  Intermodulation  Products 


A  more  severe  requirement  results  when  considering  the  three-tone  intermodulation  products 
given  by  fi  -f-  fa  -  fa  and  fa+fa-fi ,  as  shown  in  Figure  4-40.  For  this  situation,  similar  analysis  as  above  yields 
the  requirement  that  for  equal  intensify  input  frequencies: 


n,  =  3  j 


desired 


(4-75) 


For  the  example  above,  this  would  result  in  a  maximum  diffraction  efficiency  of  3%  per  tap  frequency. 
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Figure  4-39.  Two-Tone  Intermodulation  Dynamic  Ranges  of  Longitudinal 
Te02  Bragg  Cells  for  Severed  Acoustic  Power  Densities 
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SOURCE:  Reference  21 

TR-92-SI01  0029 

Figure  4-40.  Three-Tone  Intermodulation  Products 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


Many  issues  regarding  design,  analysis,  and  fabrication  have  been  presented  in  this  final  report. 
Careful  examination  of  the  adaptive  system  in  light  of  the  radar  requirements  is  necessary  to  develop  such 
issues  so  they  are  not  surprises  in  the  future.  The  synergistic  combination  of  processor  simulation  and 
hardware  development  should  offer  greater  insight  into  the  key  aspects  of  the  system.  For  example,  the 
incremental  changing  of  the  adaptive  weights  using  "block  LMS"  techniques  must  be  further  examined. 
There  are  a  number  of  challenges  yet  to  overcome,  but  we  are  confident  that  the  continued  development  of 
this  system  will  yield  elegant  solutions  that  take  full  advantage  of  the  best  attributes  of  optics  and 
electronics.  It  is  hoped  that  this  serious  consideration  of  algorithm  and  hardware  analysis  and  design  will 
direct  Dynetics’  resources  to  achieve  the  best  possible  product  in  the  long  run. 

5.1  RADAR  INSERTION  CONSIDERATIONS 

It  appears  that  the  MADOP  is  maturing  to  a  point  where  it  can  be  tested  in  a  radar  test-oeci.  To 
this  end.  there  are  several  issues  related  to  potential  performance  and  system  implementation  that  must  be 
addressed.  One  of  these  is  a  concern  for  how  the  MADOP  performs  relative  to  a  multichannel  jamming 
suppression  system  that  does  not  contain  the  delay-line  feature.  A  further  concern  is  whether  or  not  the 
existing  version  of  the  MADOP  offers  better  performance  than  an  interference  canceller  that  does  not 
include  the  delay-line  feature.  There  is  a  related  Issue  of  whether  the  nondelay-line  version  could  be  easily 
implemented  with  optical  processing  techniques.  It  is  recommended  that  this  area  be  pursued  with  analysis 
and  simulation. 

The  C-band  radar  operated  by  OCDR  appears  to  offer  the  most  attractive  Rome  Laboratory  test¬ 
bed  for  the  MADOP.  All  the  receiver  channels  are  readily  available  and  provide  outputs  at  the  proper  IF. 
They  also  appear  to  have  bandwidths  that  are  compatible  with  the  capabilities  of  the  MADOP.  The  radar 
has  a  second  auxiliary  antenna  that  could  be  mounted  near  the  existing  auxiliary  antenna  to  provide  two 
auxiliary  channels.  Furthermore,  it  might  be  possible  to  form  additional  auxiliary  channels  by  accessing 
individual  subarrays  of  the  main  array.  However,  before  such  an  approach  is  tried,  it  would  be  necessary 
to  obtain  details  on  the  scanning  characteristics  of  the  subarray.  If  the  subarrays  are  scanned,  this  could 
introduce  complications  into  the  operation  of  the  MADOP.  Other  attractive  features  of  the  C-band  array 
radar  are  the  two  jammers  and  a  boresight  tower  associated  with  it.  If  a  corner  reflector  or  some  type  of 
signal  repeater  could  be  attached  to  the  boresight  tower,  it  could  provide  a  target  signal  source  for  some  of 
the  testing. 

Another  area  of  performance  that  should  be  considered  is  the  convergence  time  of  fhe  MADOP. 
This  impacts  the  implementation,  use  and  performance  of  the  radar  when  it  incorporates  the  MADOP.  It 
may  be  that  all  pulses  transmitted  during  the  convergence  period  will  not  be  available  for  other  signal 
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processing  functions.  The  net  effect  of  this  is  a  loss  in  radar  resources.  If  the  convergence  period  is  a 
significant  portion  of  the  processing  interval,  this  could  seriously  impact  system  performance.  Another  area 
related  to  convergence  time  is  the  use  of  the  MAOOP  in  a  continuously  scanning  search  radar.  If  the 
convergence  time  is  too  large,  the  scan  rate  of  the  radar  would  have  to  be  slowed,  which  could  affect  overall 
performance  of  the  radar.  If  jamming  is  severe  enough,  the  resultant  performance  degradation  could  be 
much  worse  due  to  very  long  convergence  times.  Also  related  are  tne  effects  of  MADOP  operation  on 
clutter  cancellation,  pulse  distortion,  and  target  parameter  measurements. 

In  the  area  of  system  implementation,  the  problem  exists  as  to  just  how  the  jammer  suppression 
system  should  be  implemented  in  the  sum  and  difference  channels  of  the  radar  and  what  effect  the  jammer 
suppression  system  would  have  on  monopulse  performance.  It  should  also  be  noted  that  for  search  radars, 
which  typically  do  not  employ  monopulse  techniques,  this  may  not  be  a  problem;  for  track  radars  however, 
it  will  definitely  be  an  issue.  Furthermore,  the  impact  on  performance  will  depend  upon  the  type  of 
monopulse  processing;  e.g.,  two  channel,  three  channel,  or  four  channel.  There  is  also  a  need  to  consider 
the  implementation  impact  on  multiple-beam  antennas. 

In  Section  2,  we  discussed  how  the  MAOOP  technology  could  be  inserted  into  fielded  systems 
and/or  systems  under  development.  In  the  future,  it  will  be  necessary  to  determine  who,  currently  within  the 
Air  Force,  is  designing  new  radars  and  retrofitting  or  upgrading  older  radars.  The  people  identified  should 
be  those  who  are  in  a  position  to  authorize  or  recommend  inclusion  of  new  subsystems  in  the  radars.  A  key 
facet  to  describing  the  capabilities  is  to  define  the  MADOP  in  terms  that  are  meaningful  to  a  radar  person; 
this  includes  careful  use  of  terminology.  It  is  recommended  that  performance  results  based  on  waveforms, 
modes,  antenna  scan  characteristics,  etc.,  that  are  compatible  with  the  potential  customer’s  radar  be 
provided.  Cost  and  schedule  estimates  for  one  of  the  MADOP  systems  should  also  be  determined.  Finally, 
if  the  potential  customer  shows  any  interest,  he/she  should  be  invited  to  Rome  Laboratory  to  observe  a 
demonstration  of  the  MADOP  in  one  of  the  Rome  Laboratory  radar  test-beds  when  such  a  demonstration 
milestone  is  reached. 

5.2  MADOP  HARDWARE  DEVELOPMENT  CONSIDERATIONS 

As  summarized  in  this  final  report,  progress  on  the  design  and  hardware  implementation  of  the 
two  AO  subsystems  has  been  significant.  Both  subsystems  are  compatible  with  processing  at  the  radar 
system  IF  and  allow  for  coherent  (magnitude  and  phase)  processing.  Two  architectural  variants  of  each  AO 
subsystem  were  examined  in  both  theory  and  hardware,  and  evaluated  to  assess  ultimate  system  utility. 
The  PC  interface  between  the  two  AO  subsystems  remains  a  key  design  issue  and  needs  to  be  further 
emphasized. 
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The  testing  performed  on  the  in-line  correlator  has  given  much  insight  to  the  architecture  setup, 
operation,  capabilities  and  limitations  as  described  in  Subsection  4.2.1 .  Further  analysis/testing  should  be 
performed  since  operation  and  performance  of  the  new  components,  such  as  the  multichannel  AO  cell  and 
CCD  array,  may  differ  from  that  of  the  current  corr»^nents  in  the  preliminary  setup.  Due  to  time  and 
hardware  constraints,  full  characterization  of  the  in-i..ie  architecture  did  not  occur.  Considering  the 
architectural  comparisons  presented  in  Subsection  4.2.3,  it  was  decided  that  the  emphasis  for  hardware 
implementation  and  testing  would  be  the  two-path  correlator  architecture.  It  is  recommended  that  work  on 
the  in-line  correlator  hardware  implementation  proceed  on  an  as-needed  basis. 

The  two-path  Mach-Zehnder  correlator  development  reported  on  in  Subsection  4.2.2  has  offered 
a  second  time-integrating  correlator  architecture  that  is  judged  to  have  superior  performance  over  the  in¬ 
line  architecture.  Whereas  the  in-line  architecture  employs  a  laser  diode  and  one  AO  cell,  the  Mach- 
Zehnder  architecture  uses  two  AO  cells.  For  multichannel  operation,  either  one  or  two  multichannel  AO 
cells  can  be  employed,  although  current  testing  has  been  performed  with  one  multichannel  and  one  single¬ 
channel  AO  cell.  It  is  recommended  that  two  multichannel  AO  cells  be  employed  in  future  implementations 
to  reduced  the  optical  system  complexity  and  improve  stability.  Packaging  of  this  MADOP  subsystem  on  a 
2  by  2  ft  breadboard,  and  additional  testing  should  be  emphasized.  T ests  that  should  be  performed  include 
further  dynamic  range  characterization,  linearity,  delay  window  uniformity,  and  the  effects  of  the  post- 
detection  processing  algorithms  on  these  measurements. 

Two  AOTDL  architectures  were  also  developed  to  perform  the  filtering  operations  required  in  the 
MADOP.  The  Mach-Zehnder  interferometric  architecture,  described  in  Subsection  4.3.1,  was  an  initial 
attempt  to  obtain  coherent  operation  at  the  system  IF.  Although  heterodyne  efficiency  was  high  and  output 
was  at  the  system  IF,  it  was  recognized  theoretically  and  obsenred  experimentally  that  true  time-delay  was 
not  achievable.  In  order  to  remedy  this  situation,  a  noninterferometric  architecture  was  designed  and 
implemented. 

The  noninterferometric  AOTDL  architecture  was  designed  and  implemented  as  described  in 
Subsection  4.3.2.  This  AO  system  employs  an  electronic  reference  tone  to  improve  heterodyne  efficiency. 
The  detected  signal  output  is  at  a  reduced  center  frequency  and  is  mixed  back  up  to  the  system  IF  after 
detection.  Initial  testing  of  frequency  response  and  cancellation  performance  were  accomplished  but  further 
testing  must  be  performed  in  future  efforts.  Preliminary  measurements  made  at  the  Photonics  Center 
showed  that  true  time-delay  was  achieved  with  this  architecture,  although  later  theoretical  development 
implied  that  true  time-delay  should  not  be  observed  (see  Equation  4-69).  If  the  theory  is  verified  and  true 
time-delay  is  not  achieved,  the  electronic  reference  approach  should  be  dropped.  A  third  architecture  that 
mixes  the  undiffracted  and  diffracted  beams  (see  Equation  4-64)  should  then  be  developed.  It  is  also 
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recommended  that  the  AOTDL  subsystem  be  reduced  to  a  2  by  2  ft  package  to  improve  phase  stability, 
which  is  critical  to  noise  cancellation  performed  at  the  system  IF. 

As  noted  above,  a  key  limitation  of  the  MADOP  system  is  the  current  PC  interface.  Consideration 
of  alternative  hardware  that  offers  improved  interface  and  processing  speed  should  be  undertaken.  Also 
algorithms  for  determining  weight  functions  need  to  be  further  studied  and  coordinated  with  simulation 
activities.  Finally,  future  alternative  to  the  current  PC  interface  being  considered  is  a  combined  integrating 
photodiode/laser  diode  array.  The  integrating  photodiodes  would  replace  the  CCD  arrays  in  the  time- 
integrating  correlator  while  the  laser  diode  array  would  replace  the  AOSLM  to  provide  input  weights  to  the 
AOTDL  subsystem.  Technology  for  implementing  this  approach  needs  to  be  further  evaluated. 
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APPENDIX  A.  MATLAB  SIMULATION  ROUTINES 


This  appendix  contains  software  listings  described  in  Subsection  3.2  of  this  report.  Parameters 
employed  in  these  MATLAB  routines  are  summarized  in  Table  A-1,  and  a  flowchart  for  OJC1  is  shown  in 
Figure  A-1. 


Table  A-1.  Parameters  Employed  in  OJCl,  OJC2,  and  OJC3 


Parameter 

Definition 

d 

Received  main  channel  desired  signal 

s 

Uncorrupted  received  signal 

basenoise 

Random  Gaussian  noise  vector  (filtered  in  OJC3) 

nm 

Main  channel  jammer  signal 

SJ 

Signal  to  Jammer  ratio  in  main  channel 

na 

Auxiliary  channel  jammer  signal 

AM 

Auxiliary-to-main  channel  jamming  signal  voltage  ratio 

AMI 

AM  for  Jammer  1  in  OJC2 

AM2 

AM  for  Jammer  2  in  OJC2 

pwrorig 

Initial  jammer  power  in  d 

pwr 

Residual  jammer  power  in  e 

w 

weight  vector 

delw 

Unthresholded  weight  update  vector 

delwold 

Previous  weight  vector 

absdelw 

magnitude  of  delw  normalized  between  0  and  1 

k 

Vector  of  indices  for  which  absdelw  <  0.1 

y 

Output  estimate  of  d 

e 

Error  signal,  e  =  d  -  y 

r 

Correlation  of  e  and  na  used  to  generate  weight  update 

a 

Acceleration  parameter  [a(t)  in  Equation  3-2] 

aold 

previous  acceleration  parameter 

[bl.al] 

Poles  and  zeros  of  Butterworth  filter  in  OJC3 

[rmax,  indx] 

Maximum  value  of  r  and  its  index 
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Figure  A-1.  Flowchart  for  OJCl 
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%  optical  Jarmdng  Canceller  #1 
%  "white”  noise  case 

%  Iterative  but  reuses  sane  signal  and  noise  -  not  real  world 
%  Uses  Welstead's  ad^ive  method  for  ccnpjting  a  in  w=wfa*r 
%  Only  one  jammer  and  one  multipath 
%  Uses  thresholded  weight  update 
a^ .  3; 

delfcrold=zeros (1, 100) ; 

% 

%Spacing  between  time  saitples  is  0.05  ps  (20  Mhz  sanple  rate) 
%Signal,  s(:)  is  a  2  ps  pulse  starting  at  1  ps 
%Total  duration  of  data  is  10  (is 
s=zeros(l,  400)  ; 
s  (121: 160)  =ones  (1,40)  ; 

%  Generate  jammer 

%  For  now  jammer  is  noise  with  a  BW  of  l/(0.05  (is)  =  20  MHz 

%  It  is  also  Gaussian 

rand( 'soed',0) 

rand  ('normal') 

basenoise=rand ( 1, 500) ; 

%  basenoise  is  the  basic  jammer  source  -  from  which  all  others  are 
%  derived.  It  represents  the  actual  jammer  output  (sort  of) 
%Generate  main  channel  jamming  signal  -  delay  basenoise  ty  2  |is 
nmpbasenoise (21 : 420) ; 

%  Generate  main  channel  signal  plus  jamming,  d(t) 

%  SJ  is  the  signal-to-jam  ratio  (volt/volt) 

SJ=sqrt(0.1);  %  -10  dB 

%  Corpute  jarrmer  power 

pwrorig=10*logl0(  (rin(101:300)  *nm(101:300)  •)  /200/SJ/S7)  ; 

<*=  s  +  nm/SJ; 

%  Generate  aux  channel  jamming  signal  -  delay  Uisaioise  by  1  (is 
%  AM  is  the  aux-to-rrain  ratio  (volt/volt) 

AM=sqrt(0.5);  %  (-3  dB) 

AM=1; 

na=t)asenoise  (41:440)  *AM; 

%  DefLne  initial  weights  as  all  zeros 
w=2eros  (1, 100)  ; 

%%%%%%%%%%%%%%%%%%% 

%  Start  loop 
%%%%%%%%%%%%%%%%%% 
for  iter  =  1:100 

%  form  y=Sw*na  -  the  weighted  aux  channel  signal 
y=f liter  (w,  1, na)  ; 

%form  error  signal,  e=d-y 
e=d-y; 

%  corpute  output  power  after  cancellation 
res=e ( 101 : 300) -s (101 : 300) ; 
pwr=10*logl0 ( (res*res') /200) ; 

%  Correlate  e  and  na  call  the  result  r 
r=xoorr (e, na, 'biased' ) ; 

%plot  interim  results  -  error  and  correlation 
clg 

sulnplot  (221)  ,plot  (e  (101 -.300) ) ;  title  ('error') 

subplot (222) , plot (r (351 :550) ) ;title( 'correlator  output') 

%update  weight  vector 
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%  use  zesidues  above  threshold 
del*^r(400:499) ; 

afaedelwabs  (delw)  /max  (abs  (delw) )  ; 
lc->find(absdelw<0.1) ; 
delw  Oc)  -zeros  (k)  ; 
numer*<jelwold  *  delwold' ; 
donanKielwold  *  (delwold-delw) ' ; 
if  (nunner*denan)  ■=  0 
a=0.5; 
else 

a=aold*nLimer/ciaxnv 

end 

aold=a; 

delwold=clelw; 

v«ri-a*delw; 

%  pick  peak  weigfit  only 
%tnTBx, indx]=iTBx(r(400;499) )  ; 
%delw=zeros  (1, 100) ; 

%del  w  ( indx)  =rmax; 

%nuter=delwold  *  delwold'; 
%dencin=delwold  *  (delwold-delw) '  ; 
%if  (numer*dencm)  —  0 
%  a=0.5; 

%else 

%  a=aold*numer/denati; 

%end 

%aold=a; 

%delwold=delw; 

%w=w+a*delw; 


— >> 

I 

\ 

^  Threshold  Weight  Vector  Update 


%  use  all  residues 
%deiw=r(400:499); 

%nurer=clBlwold  *  delwold'; 

%denotm<jelwold  *  (delwold-delw) '  ; 

%if  (numer*daian)  =0  I 

%  a=0.5;  )  Full  Weight  Vector  Update 

%else  f 

%  a=aold*nuter/dencm; 

%end 

%aold=a; 

%delwold=delw; 

%w=vH-a*delw;  ^ 

subplot  (223)  ,plot  (delw)  ;title  ( 'thresholded  weicfit  increments ' ) 
subplot  (224 ) ,  plot  (w) ;  title  ( '  weic^it ' ) 
label-  ^rintf (' Iteration  =  %3.0f,  a=  %6.3f ',iter,a) ; 
labell-  sprintf  Ciip  jmr  pwr  =  %5.1f,  out  jmr  pwr-  %5.1f  (in 
dB)  ',pwrorig,pwr)  ; 
text (0.1, 0.5, label, '  C) 
text (0.5, 0.5, labell,  sc') 


end 
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%  Cptical  JannvLng  Canceller  #2 
%  "wtiite"  noise  case 

%  Iterative  txit  reuses  sarte  signal  and  noise  -  not  real  world 
%  Uses  Welstead’s  adaptive  method  for  carputing  a  in  w=vri-a*r 
%  Only  cne  jartiter  and  but  two  nultipaths  in  each  antenna 
%  Uses  thresholded  weight  ipdate 
a=.3; 

clslwDld=zeros  (1, 100)  ; 

% 

%Spacing  between  tirte  sanples  is  0.05  ps  <20  Mhz  sanple  rate) 
%Signal,  s  (t)  is  a  2  )is  pulse  starting  at  1  ps 
%Total  duration  of  data  is  10  ps 
s=zeros  (1, 400)  ; 
s  (121: 160)  =cxies  (1,40); 

%  Generate  jairmer 

%  For  now  jammer  is  noise  with  a  EW  of  l/(0.05  ps)  =  20  MHz 

%  It  is  also  Gaussian 

rand  (' seed ',0) 

rand  ('normal') 

basenoise=rand (1,500) ; 

%  basenoise  is  the  basic  jamter  source  -  from  v^ch  all  others  are 
%  derived.  It  represents  the  actual  jammer  output  (sort  of) 
%Generate  main  channel  jamming  signal  -  delay  basenoise  by  2  ps 
nm=basenoise(21:420)  +  0.5*basenoise(31:430) ; 

%  (Generate  main  channel  signal  plus  jamming,  d(t) 

%  SJ  is  the  signal-to-jam  ratio  (volt/volt) 

SJ=sqrt(0.1) ;  %  -10  dB 

d=  s  +  nm/SJ; 

%  (jenerate  aux  channel  jamming  signal  -  delay  basenoise  by  1  ps 
%  is  the  aux-to-main  ratio  (volt/volt) 
flMl=sqrt<0.5) ;  %  (-3  dB) 

AM2=1;  %  (OdB) 

na=toasenoise(41:440)*Wtl  +  basenoise  (66: 465) 

%  Define  initial  weights  as  all  zeros 
w=2eros (1, 100) ; 

%%%%%%%%%%%%%%%%%%% 

%  Start  loop 
%%%%%%%%%%%%%%%%%% 
for  iter  =  1:100 

%  form  y=5>r*na  -  the  weighted  aux  channel  .signal 
y=f liter (w, 1, na) ; 

%form  error  signal,  e=d-y 
e=d-y; 

%  Correlate  e  and  na  call  the  result  r 
r=xoorr (e, na,  'biased') ; 

%plot  interim  results  -  error  and  oorrelaticn 
clg 

subplot (221) , plot (e (101:300) ) ;title ( 'error' ) 

subplot (222) ,plot (r (351 :550) ) ;title ( 'correlator  output' ) 

%update  weight  vector 
delv«^r  (400:499)  ; 

absdelw=abs  (delw)  /max.  (abs  (delw) )  ; 

)c=find(absdelw<0.1)  ; 
delw  (k)  =zeros  (k)  ; 
numer=clelwold  *  delwold'; 
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clenan*<lelwold  *  (delwoldhclelw) ' ; 
if  (mner^denan)  —  0 
a-0.5; 
else 

a*^old*nuier/<jaKin; 

end 

aold=a; 

CElwoldKtelw; 

»«#+a*dBlw; 

subplot (223), plot  (delw)  .-title Cthresholded  weight  increments') 

subplot  (224)  ,plot  (w)  .-title  ( 'weight ' ) 

label=  ^rintf (' Iteration  =  %3.0f,  a=  %6.3f',iter,a).- 

text (0.45, 0.5, label, 'sc') 

end 


%  Optical  Janndng  Canceller  #3 

%  bandlimlted  noise  case  -  noise  filtered  by  2  Miz,  5  pole  Butterwortb 
%  filter.  Iterative  but  reuses  same  signal  and  noise  -  not  real  world 
%  Uses  Weistead's  ad^ive  method  for  oonputing  a  in  w=^a*r 
%  Cnly  one  jammer  and  one  nultipath 
%  Uses  thresholded  vgeig^t  vqpdate 
a=.3; 

delwold-zeros (1, 100) ; 

%0pacing  betweai  time  sanples  is  0.05  ms  (20  Mhz  sanple  rate) 

%Slgnal,  s(t)  is  a  2  ms  pulse  starting  at  1  ms 
%Total  (iiration  of  data  is  10  ms 
s>°2eros(l,  400)  ; 
s  (121: 160)  =ones  (1,40); 

%  (Generate  jammer 

%  For  now  jammer  is  noise  with  a  BW  of  l/(0.05  Ms)  =  20  MHz 

%  It  is  also  (Russian 

rand  (' seed',  0) 

rand('nocnal') 

basenoise-rand  (1,500) ; 

[bl,al]Hxitter(5,2/10) ; 
basaioise-filtfilt (bl,al,basenoise) ; 

%  basenoise  is  the  basic  jammer  souroe  -  from  which  all  others  are 
%  derived.  It  represents  the  actual  jammer  output  (sort  of) 

%Generate  main  channel  jamming  signal  -  delay  basenoise  by  2  ms 
nn43as»K)ise  (21 : 420)  ; 

%  Generate  main  channel  signal  plus  jamming,  d(t) 

%  SJ  is  the  signal-to-jam  ratio  (volt/volt) 

SJ-9qrt(0.1);  %  -10  dB 

%  (3cnpute  jammer  power 

pwrorig-lonoglO  ( (nmdOl  :300)  *nm(101 :300)  • )  /200/SJ/SJ)  ; 
s  +  nm/SJ; 

%  Generate  aux  channel  jamming  signal  -  delay  basenoise  by  1  Ms 
%  Ml  is  the  aux-to-main  ratio  (volt/volt) 

Ml-sqrt(0.5);  %  (-3  dB) 

M4-1; 

na^sasenoise  ( 41 : 440)  *M1; 

%  Define  initial  weights  as  all  zeros 
f^>zetos  (1, 1(X)) ; 

%%%%%%%%%%%%%%%%%%% 

%  Start  loop 
%%%%%%%%%%%%%%%%%% 
for  iter  -  1:100 

%  form  y-£w*na  -  the  weighted  aux  channel  signal 
y*f  liter  (w,  1,  na) ; 

%form  error  signal,  eri-y 
e^y; 

%  Ccnpute  output  power  after  cancellaticn 
tes-«  ( 101 : 300)  -s  (101 : 300)  ; 
pwr-10*logl0  ( (res*res ' )  /200) ; 

%  Correlate  e  and  na  call  the  result  r 
r>TOorr  (e,na,  'biased' )  ; 

%plot  interim  results  -  error  and  correlation 
rig 

svhplot  (221) ,  plot  (e  (101 : 300) )  ;title  ( '  error' ) 

subplot (222), plot (r (351:550))  /title ('correlator  output') 
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%update  weight  vector 
%use  residues  above  threshold 
%delw=r (400:499); 

%absdelw=ebs  (delw)  /max  (abs  (delw) )  ; 

%k=f  ind  (absclelw<0 . 1 ) ; 

%delw  (k)  “zeros  (k)  ; 

%nurer=<telwold  *  delwold'; 

%daicm=delwold  *  (delwold-delw) '  ; 

%if  (nLStBr*cteiam)  =  0 
%  a“0.5; 

%else 

%  a»aold*nuter/dencnv; 

%end 

%aold=!a; 

%delwold=<ielw; 

%v«irt-a*delw; 

%  pick  peak  weight  only 
%[nnax,indx]*TOx(r(400:499) )  ; 

%delvp=zeros  (1, 100)  ; 

%delw  (indx)  “max; 

%numer=delwDld  *  delwold'; 

%denonK]elwold  *  (delwold-delw) '  ; 

%if  (numer*denan)  =  0 
%  a“0.5; 

%else 

%  a“aold*niitBr/dencitu 

%end 

%aold=a; 

%delwold-delw; 

%v«irt-a*delw; 

%  use  all  residues 
delw»r(400:499) ; 
numBr=delwold  *  delwold'; 
denanpdelwold  *  (delwold-delw) '  ; 
if  (nuner^denati)  “=  0 
a=0.5; 
else 

a“aold*ntBtBr/denom; 

end 

aold>«; 

delwold-delw; 

w=^a*delw; 

%  use 

svtoplot  (223)  ,plot  (delw)  ;title  ( 'thresholded  weight  increnents ' ) 
sxiplot  (224)  ,plot  (w)  ;title  ( 'weight ' ) 
label”  ^rlntf (' Iteration  “  %3.0f,  a  =  %6.3f ',iter,a); 
labell”  sprintfCinp  jmr  pwr  “  %5.1f,  out  jmr  pwr=  %5.1f  (in 
dB)  ',pwrorig,pwr)  ; 
text (0.1,0. 5, label, 'sc') 
text (0.5, 0.5, labell,  'sc') 

end 
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MISSION 


ROME  LABORATORY 


Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search,  development,  test,  and  technology  transition  in  support  of  Air 
Force  Command,  Control,  Communications  and  Intelligence  (C^I)  activities 
far  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  "within 
areas  of  competence  is  provided  to  BSD  Program  Offices  (POs)  and  other 
BSD  elements  to  perform  effective  acqui^tion  of  C^I  systems.  In  addition, 
^dime  Labaratorys  technology  supports  other  AFSC  Product  Divisions,  the 
Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  c<^petence  arui  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  inf(^mation  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainability  and  testability. 


